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PREFACE 


This report presents the results of studies conducted during 
the period July 1, 1969 - October 31, 1970, under NASA research contract 
NAS 8-25102, "Mathematical Characterization of Mechanical Behavior of 
Porous Friction Granular Media". This study was monitored by 
Dr. N. C. Costes, the Geotechnical Laboratory of NASA's Marshall Space 
Flight Center. 

The objectives of this project are: 

(1) to develop a consistent three-dimensional mathematical theory de- 
scribing the mechanical behavior of porous, frictional granular media 
exhibiting a small amount of cohension. — Volume I of this report 

by Dr’. G. Aquirre-Rameriz . 

(2) to solve boundary-value problems related to in-situ" measurements 
performed on the lunar or planetary surface. -- 
report by Dr. T. J. Chung and Mr. J. K. Lee; 
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SECTION I 

ON THE SOILD MASS CONTINUUM 
1. INTRODUCTION 

A general accepted practice for the analysis of the response of a 
soil mass under external loads is to model the soil mass as a continuum. 

Under this assumption the distinction between modeling of a soil mass and 
a strictly solid mass (such as rolled steel, say) as a continuum is made 
through the constitutive equations. 

In this paper the basic equations of continuum mechanics are reassessed 
as to their . applicability to model a soil mass as a continuum. For the 
purpose of simplicity the discussion is limited to a dry (i.e., unsaturated) 
soil mass. Under these conditions the effect of interaction in a water 
soil system need not be considered. . . 

As a starting point the two physical properties that bodies are known 
to have and which are used in the construction of a model for continua-are 
taken as fundamental. These are Cl] (a) that they occupy regions of space, 
and (b) that they have mass. These two porperties are used together with the 
porous geometry' of a soil mass to show that in order to define the thermo- 
dynamic state variable *p, the density of solid aggregate, as a field variable 
one needs to introduce two field variable p(x,t), the soil bulk mass density, 
and n(x,t) the porosity of the soil. By interpreting the mass density appear- 
ing in the local balance laws of mass, momentum, and energy of continuum 
mechanics as the bulk mass density, these balance equations can be taken over 
unchanged into soil mechanics to locally describe the corresponding balance. 


laws for the soil. 
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In solid mechanics the two soil geometric variables, surface porosity 
■and volume porosity are taken to be the same. Here it is. shown that by 
defining Terzaghi's effective stress [5] by boundary conditions this 
assumption can be removed. Researchers in soil mechanics have consistently 
verified that the mechanical behavior of the effective soil structure is 
governed by the effective stress (cf. Schofield and Wroth [ll]). For this 
reason the balance of momentum and energy are formulated in terms of e-ffec- 
tive stress. 

The introduction of the porosity as a field variable requires an 
additional constitutive equation for the rate of change of porosity. The 
need for- this constitutive equation is shown to be given by the Second Law 
of Thermodynamics. It is shown that porosity influences the free energy 

of the soil as an internal state variable and since its rate of change does 

/ 

not appear in the Clausius-Duhem inequality a constitutive equation for it 
must be postulated. As an example a possible set of constitutive equations 
is examined. The restrictions on the proposed constitutive response 
functions are also found by a method introduced by Coleman and Gurtin [6] 
constitutive equations are further linearized.. 

NOTATION 

In this paper direct tensor notation is used. Second order tensors and 

linear transformations of the three dimensional vector space U into itself 

T 

are regarded as the same. If T is a linear transformation, T denotes its 
transpose, T _1 its inverse, tr T its trace and det T its determinann. 

The gradient with respect to spatial coordinates is denoted by grad and 
the gradient with respect to material coordinates by V. 
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2. THE SOIL MASS CONTINUUM 

It is important to realize at the outset that the' two. physical 
properties that bodies are known to have and which are taken in the 
construction of a model for continua are Cl]: (a) that they occupy regions 
of space, and (b) that they have mass. These two properties will now be 
examined in relation to a soil mass. 

Let R be an arbitrary region of space and consider two bodies B^ and 
B 2 with masses and respectively. B^ will be considered to be a 
strictly solid body (rolled steel, say) but B 2 will definitely be taken as 
a dry soil. Assuming for the present that B^ and 3 ^ occupy equal regions of 
space R, densities p^ and can be introduced such that 



R 


( 2 . 1 ) 



‘ (2.2) 


The density p^ for so introduced is the mass density . However, since B^ 
has mass by virtue of its soil skeleton, the density p^ is the bulk mass 
density . Thus, for the solid continuum and soil mass continuum the defini- 
tions (2.1) and (2.2) yield two different types of mass densities. 

The question arises as to whether a mass density for the soil body can 
be introduced as a field variable. The answer to this question is positive- 
provided that another field variable be admitted. To this purpose let V g 
denote the' volume of the solid aggregate of the soil body 3 ^ which occupies 
the region of space R. A density UJ can then be introduced such that 



(2.3) 
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The density W is the additional field variable that may be .introduced for 
the soil mass continuum [ 2 ] which, if it is also introduced for the solid 
continuum, has a constant value of one. The variable ui may be called the 
solids material volume density. Note that if dv is an element of volume 
of soil mass then 

dv^ ='U)dv .(2.4) 

gives the element of solid aggregate volume. Having introduced U) the more 
familiar variable, porosity,- of the soil mass continuum is obtained from 

n = 1 - U). (2.5) 

The mass density p for the soil mass continuum is then given by 



Therefore, the soil mass continuum differs from the solid continuum in 
the sense that two field variables, p^, W are needed to obtain the mass' 
density of the solid aggregates. This idea will now be put into formal 
grounds. 

'In continuum mechanics a body B is considered to be a manifold of 
particles, denoted by X. The particles, however, are primitive elements 
in the sense that numbers are primitive elements in analysis [ 3 ]. The 
-body manifold is further assumed to be (1) smooth and isomorphic to regions 
in Euclidean 3-space, and (2) endowed with a non-negative measure M of space 
which is its mass distribution [3]. In the same, sense, a soil body can be 
considered to be a set of pnrtic]os, denoted by X. These particles may be 
called 'soil particles which, of course, are not to be confused with the 
physical soil particles. The soil body manifold can also be assumed to be 
(1) smooth and isomorphic' to regions in Euclidean 3-space, and (?) endowed 
with two non-negative- measures of space: (a) M, which is the distribution 
of mass of the .solid aggregate, and (b) V g , which is the distribution of 
material volume of the solid aggregate. The measure M is assigned once . 



and for all. This is not so with the measure V which can vary in time 

s 

due to local microscopic deformation of the solid aggreagre. 

Let B denote a soil body. B can then be considered a set of particles 
X, with the above mentioned structure. A configuration of B is then a 
smooth homeomorphism ^ of B onto a subset of Euclidean 3-space £: 

x = X(X), X = X _1 ( x ) • (2.7) 

Here X ^ denotes the inverse of the mapping X* 

Consider a small element of cross-sectional area A A of the soil body 
in the configuration To an observer AA will then appear -to have a 
"Swiss cheese- like" structure in which part of AA will correspond to the 
solid aggregates and the remaining part will be voids. Assuming that the 
voids are randomly distributed in their location within AA then in the limit 
as AA becomes infinitesimal a distribution CL (x) can be assigned to the 
ratio of dA v> the element of area occupied by the voids, to dA, i.e., 

dA v =^(£>dA. (2.8) 

Equation (2.8) will be used when the soil stresses are discussed. 


3 . KINEMATICS 



A motion of the soil body B is a one-parameter family of configurations 


x 

rw 


£ t (x) - xcx.t ). 


(3.1) 


where the parameter te(-“, co ) is the time. In view of -the smoothness 
assumptions made )( will be invertible for each te(- co , 03 ). 

A reference configuration for B is a fixed configuration^. The place 
of X £ B in its reference configuration K is denoted by X, 
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Insertion of (3.2)^ into (S.l)^ leads to 

£ = &<£>'> -tf* 1 ©- c >- < 3 - 3 > 

■ ~ t 

The function X is the deformation function for B from its reference 

fvK 

configuration M.. 

The velocity and acceleration for XEB are defined by 


x 


a_ 

at 




X 

rv 



X R (X,t) 

rwH rw 


while the deformation gradient is given by 

£- v V£. c ) 

/V 


(3.4) 


(3.5) 


F is a second order non-singular tensor with the property 

|det f| > 0. (3.6) 

If L denotes the gradient of the velocity x then it can be shown that 

rs/ 


■ L 



= grad x(x,t) . 

'v fs/ 


(3.7) 


The various strain measures used in continuum mechanics are constructed 
from the deformation gradient F [4] and these can be used to describe the 

rv 

deformation of the soil mass continuum. However, for the soil mass continuum 
an additional "strain" measure can be constructed'. If dv denotes an 
element of volume' in ft(B) at X and dv its image under the mapping (3.3) 
then it is known that [ 4 ] 

dv = |det F|dV\ (3.8) 

Letting dV be the clement of solid aggregate volume in R(B) and dv the 
s ^ s 

corresponding element in^ ?t (B,t) the following quantity can be constructed 
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dv (bidet f[ 
s _ 1 ~ 

dv ” UJ 
s o 


(3.9) 


where W is the solid aggregate volume density in w.(B). A is the average 

O ^ 

expansion (or contraction) of the solid aggregate. The observation is made 
that if the assumption of incompressibility of the solid aggregate of the 
microscopic level is introduced, as is common practice in soil mechanics, 
then 

A = 1 . (3.10) 


and (3.9) becomes 


U) = wldet f|. 

o 1 ~ 


(3.11) 


Introducing the change in porosity, a through 

a, = n - n 

O 


(3.12) 


where n Q is the porosity in H(B), Eq. (3.9) .can be written in the form 

A = (1 - y-™~) jdet fJ . (3.13) 

o 

Whenever the deformation of the soil mass continuum is infinitesimal the 
linear strain measure E given by 

E = -^(H + h T ) , . (3.14) 

whore ]l i s Lhe displacement gradient, can be used. In this case under the 
assumption of small displacement gradients, i.c., /trini^ « 1, 

jdet F | ^ 1 + tr E. ■ (3.15) 


Then assuming Ot is of the same order as the strains, E, Eq. (3.9) can be 
written as 


A - 


1 = 


tr IT 


1 


a 


- n 

o 


(3.16) 
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Tlic term tr ^ is the volumetric strain. When the solid aggregate is assumed 
to be microscopically incompressible (3.16) becomes 


tr E 



(3.17) 


which is the formula used in soil mechanics to compute infinitesimal 
volumetric strains. 

Quite often, in soil mechanics, the void ratio e is used instead of 
the porosity. This is defined by 


e 




(3.18) 


where civ is the. element of volume of the voids. Note that if e and e 
v o 

denote initial and current void ratios respectively then 


Q- = C 

1 - n 1 + e 
o o 


(3.19) 


where £ = e - e . Equation (3.17) can then be. written in the form 


tr IT 


__C 

1 + e * 
o 


(3.20) 


4. THE BALANCE EQUATIONS 

The .basic balance equations of continuum mechanics arc local formu- 
lations of the. principles of physics of conservation of (a) mass, (b) linear 
momentmn, (c) moment of momentum, (d) energy. These will be taken one at a 
t ime . 

Balance of Mass - The balance of mass equation is given in one of two 
forms [ 4 ]: the. spatial form 


-I 7 '+ div p x = 0 


(4.1) 
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or the material form 


P |dct fJ = P q (4.2) 

where p and p are the mass densities in the reference and current configu- 
o 

rations respectively. 

It can be shown that (4.1) and (4.2) can be taken over unchanged as 
the balance equation for the soil mass continuum provided that the mass 
density appearing therein be interpreted as the "bulk mass density". 
However, for the soil mass continuum (3.8) is available. This equation 
can be .used to eliminate jdet f} from (4.2) and obtain 

p A = P o (4.3) 


for the balance of mass equation in terms of "mass density" of the solid 
aggregate. 

Balance of Momentum and Moment of Momentum - Consider an element of oriented 
surface area da, of the soil mass continuum in its current configuration and 
denote the area fraction defined by (2.8) by 

a = a(x,t) = d . ■ (x) (4.4) 

4>t 

Further lot t be the stress vector acting on da. The stress vector t 
— ~ ~s 

acting on the solid aggregate- portion of d£. is defined by 

tg = at. (4.5) 


Cauchy's stress hypothesis is invoked and the existence of a stress tensor 
T is assumed such that 

<-vS 


t = T n 


(4.6) 


where n is the unit normal vector of orientation of dOk. The stress tensor 
T is called Terzaghi's effective stress. The reason for calling T the 
effective stress will now be given. 
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Consider a completely saturated soil and let T, be the total stress 
acting at a point of the soil'. Terzaghi's effective stress is then 
defined by [ 5 ] 

T = T + pi (4.7) 

A/ 

where p is the pore water pressure which here is taken as positive in 
tension. Now the portion d d c of the oriented element of surface d Cu 
occupied by. the fluid is given by 

d £f = (i-^da (4.8) 

and the portion of the surface vector t acting on the fluid is given by 

' r*j 

t, = (l-0)t. • (4.9) 


Note that 


t + t c = t. 


(4.10) 


Cauchy's stress hypothesis is invoked and the existence of a total stress 

tensor X and partial stress tensor T, is assumed such that 
~ 1 ~f 


t = Tn, 


= T,n. 

~f~ 


(4.11) 


•Substitution of (4.6), (4.11) into (4.10) leads to 


(1 


■~s 


+ T 


- T) n = 0. 


(4.12) 


If (4.12) is to hold for arbitrary n it follows that 


T - T 


+ V 


(4.13) 


Equation (4.7) follows from (4 .13) by setting T^. .= pi. Thus the reason for 
calling T the effective stress. 
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Having defined the soil stress, the balance of momentum and moment of 
momentum for the soil mass continuum is then postulated to be given by 


Px = div T + pb 

rss rvS 


T = T 

~s • ~sr 


(4.14) 

(4.15) 


.wherein b is' the specific body force density, i.e., the body force per 
unit solid aggregate, mass. 

Balance of Energy - Let £ denote the specific internal energy of the soil 

body, i.e., the internal energy per unit solid aggregate mass, r the energy 

source due to external radiation, and q the heat flux. In a manner similar 

to that used, to define the effective stress, an effective heat flux q can 

~s 

be defined by 

As = 4a- < 4a6 > 


The balance of energy for the soil mass continuum is then postulated to be 
given by 

pe = tr T L - div q -i- pr. (4.17) 


5. THE CLAUSIUS- DU11KM INEQUALITY 

Let 0 = ®(x,t) be the temperature of the soil mass continuum which is 
assumed to be positive and let T| be the specific entropy. Then regarding 
q /0 as the flux of entropy due to heat flow and r/0 the supply of entropy 
from radiation, the specific rate Y of entropy production for the soil 
mass continuum is postulated to be given by [ 4 ] 

PY = P*1 - (*§ - div 5 q s ). . (5.1) 

The Clausius-Duhom inequality is the assertion that the rate of entropy 
production is not negative, i.e.. 


Y a 0. 


(5.2) 



12 


Equation (5.1) may be combined with (4.17) so as to obtain 

. <5.3) 

pt) 

where (3.7) has been used and 

■£ “ srad 6 - L ‘ j istf)' 1 (5.4) 

The specific free energy may be introduced through the definition 

♦ “ e - 0*n. . ' - (5.5)' 

Under this definition ‘(5.4) may be written' in the form 

0Y “ - 110 + tr S F - ~ q • q. (5.6) 

The Clausius-Duhem inequality (5.2) will be used .to find restrictions on 

constitutive equations. 


6. CONSTITUTIVE EQUATIONS 

It’ is well known that the deformation of a soil mass continuum is 
accompanied by dissipative effects which are in addition to heat conduction. 
Therefore the constitutive equations for the soil body must be such as to 
show this feature. In continuum mechanics there are various ways of 
accounting for dissipative effects which (in addition to heat conduction) 
accompany deformation. One of- these is to postulate the existence of 
internal state variables which influence the free energy and whose rate of 
'change is governed by differential equations in which the strain appears. 
These have been studied in detail by Coleman and Gurtin [6] for single 
continua and by Bowen t-7,8] for mixtures of continue. It will be shown 
below that by considering the porosity of the. soil mass' continuum as an 
internal state variable the dissipative, effects can be accounted for. 
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It is important to realize at the outset that a thermodynamic state 
variable for a soil mass continuum is the solid aggrega tc mass density and 
not the bulk mass density . The solid aggregate mass density is given by p. 


The balance laws suggest that constitutive equations are needed for 

A Tl T and q . The soil mass behavior which will be studied here is one 

9 19 ~ s ~s _ _ _ 

which is characterized by four response functions T|, £, jq which give 

tlr, ri, T , and q when p, F, F, 0, g are known, i.c.. 


(f, ^ > sJ = f (P> F > L> 0 > s) 


( 6 . 1 ) 


where f stands for any q. The reason for including F as an 

independent variable will be explained below. 


Now through (2.5), (4.2), equation (2.6) may be written in the form 


P = 


(1-n) |det fJ' 


( 6 . 2 ) 


In view of this equation the existence of functions f such that 

(M, T . q ) * f(F, F, 9, g, n) (6.3) 

can be reasoned. This demonstrates how the porosity enters as an independent 
variable. Since the rate of change n of porosity does not appear in the 
Clausius-Duhem inequality a constitutive equation for n must be postulated. 
Thus to (6.3) must be added 

n = ^(F, F, 9, g, n) . (6.4) 


Therefore the porosity influences the free energy and its rate of change is 
postulated by a differential equation in which the strain appears through F. 
Of course the assumption has to be made that n, F, F, 0, g as functions of 
Xand t are smooth enough to insure the existence of a unique solution 


n = n(X,t) of (6.4) for all t in some interval [t , t + T ] with n(X, t ) = 


n (X). 
o ^ 


The inclusion of F as an independent variable will now be explained. 
It can be shown that whenever the solid aggregate is incompressible the 
porosity is governed by the differential equation 
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n = (1-n) tr F F~ 1 (6.5) 

Equation (6.5) is a special case of (6.4) with 

A A 

dnl) = 0, Bgn = 0 

U ** tsi 


so that 


'nCF, Fj n) = (1-n) tr F F _1 . 

'/v »v /V Ai 


Thus the reason for including F as an independent variable. 

It is also assumed that the effective stress T is the sum of a . 
non-dissipa tive part T^ and a dissipative part T^, i.e.. 


T - T 

As/S 


+ Id 


such that 


T = T (F, 6, n) 

fUO rv 


lo *&>&£> e >£- n > 


Therefore 


T (£» 9, *0 “£,<£»■ 9 . n > + ^(F,.?, 0, g, n) 

A/ rw rv ^ rv /V 


( 6 . 6 ) 


(6.7) 


( 6 . 8 ) 


In Section 7 it will be shown that 

A • . ' . 

i’ (F, 0, 9, 0, n) = 0 whenever n = 0. (6.9) 

/Vjj A* 

A A 

The constitutive response functions f and n have to satisfy the Clausius 
Duhem inequality (5.6). The restrictions on the response functions £ and 
by the Clausius-Duhem inequality can be found by the method used by Coleman 

and Gurtin [63 and Bowen C7,8]. The main results of this exercise are: 

A A 

(I) The response functions and r| are independent of F_ and g, i.e., 

'(^ TO « h(F, 0, n) 


(6.30) 
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A A y\ 

where h stands for either l|f or rj. 

, A A 

(II) v determines T] through the entropy relation- 

n = Q, n). (6.11) 

A A 

(III) V determines T_ through a stress relation 

^4/ 


A 


£, " <V& e > ">£ • 


(6.12) 


\\ /S Z' 

(IV) V, n, T^, and obey the general dissipation inequality 


tr &<£. i 0 - s* ">i - a „*(£. 0 > *)*<& i «. & »> 
- i- 0 > &> ■>> • s * °- 


(6.13) 


where 


Sd< 


1 A 

) = — T ( 
' p M)' 


T -1 

)<F A ) \ 


(6.14) 


Equations (6.10) through (6.13) are necessary and sufficient conditions 
that the Clausius-Duhom inequality be satisfied by the constitutive response 
functions. It is possible to extract additional information from the general 
dissipation inequality (6.13). This information will bq examined in the 
next section. 


7. EQUILIBRIUM STATES 

The additional information that can be extracted from the general 
dissipation inequality (6.13) is obtained for certain values of the indepen- 
dent variables. It has been found convenient to name equilibrium the state 
in which these values occur. 


The general dissipation inequality (6.13) implies that when £ = 0, 
the mechanical dissipation inequality 


“ tr ^D ( I s i> 0} £ j n) £ + V& 0> i> e > £» n > * 0 


(7.1) 


A A 


holds, and when (JF, ^ n) = (J3, 0) the heat conduction inequality 


&Z> £> 8 > tx) • g * 0 


(7.2) 
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holds. Also when (F, g) = (0, 0) the Internal dissipation inequality 


9 KF, 9, n)n(F, 0, 9, 0, n) £ 0 


(7.3) 


holds. 


* * 

It is convenient to call a triplet (F , 0 , n ) with 


$(P*. 0, 0*, 0, n*) = 0 


(7.4) 


an equilibrium state for the soil material point X. Note that if 
i(F, F } 9, g, n) = tr s' (F , F, 0, g, n)F - M(F, 0, n) n (F, F, 0, g, n) 


“9 5(F. j, 0, g, n) * g 


^ rv 


(7.5) 


then (6.13) can be written as 


i(F, F, 0, £, n) * 0. 


(7.6) 


Clearly 

0, 0*. 0,0 - 0 (7.7) 

Therefore .as a function of (F, F, 0, g, n) , Jl is a minimum at the equilibrium 
* A * * 

state (F , 0 , n ). Consequently 

d i*(F* + U, Xb, 9* + Xa, Xa, n* 4- Xd) L _ «= 0 (7.8) 

aA 


for all scalars a. d, all vectors a, and all second order tensors A, B in 
the domain of 

^ ^ 

For a function G( ) evaluated at (F , 0, 0 , 0, n ) the following 
notation is used 


+ , * 
G = G(F , 


9 , 6 


* 

0 n ) . 
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Equations (7.5) and (7.8) imply that 

tr(fs? - 3 i|i + 3 ♦ 'n + )B - 3 {tr 3 > n + ' A + 3,(n + a + 

V ~D n F n F ■ ~ 0 n 


j -s (1",+ -i A+ . 1 A+ s n 

-(o V 3gn 4- — -~T q ) - a = 0 
n _+ 0 ,> ~ ~ 

P d 


(7.9) 


From the second term in this equation and the arbitrariness of a, and d 
the following is concluded: either 


~i A+ _ -i A+ n A+ A 

3 n = 0 , 3 Q n = 0 , 3 n =0 

F ~ t) n 


(7.10) 


or 


A, 

3 r «* o. 

n 


(7.11) 


For obvious reasons (7.10) must be discarded. Therefore in view of (7.11) 
and the arbitrariness of B and a the following additional information is 
extracted from (7.9) 


A - * .* * 

S^(F , 0, 6 , 0, n ) = 0 




^(F*, 0, 6“, 0, n") = 0 


(7.12) 

(7.13) 


Thus at the equilibrium state the dissipative stress and the heat flux 
vanish. Also (7.11) reads 

^(£,0*,/)=°. (7-14) 

Equation (7.14) is called the equation of internal equilibrium? 


Coleman and Gurfin Cb] have derived an equation identical to (7.14) in 
their study of constitutive equations for which the independent variables 
are (F, 9, g, Ct) , a being an internal state N-vector. 
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It. la assumed that corresponding to each strain- temperature pair (F , 0 ) 

•k ~ 

there is exactly one porosity n such that (7.4) holds. This correspondence 

is given by the function 


* 

n 


* * 
a Q. > 


6*) 


(7.15) 


which is called the equilibrium response function for n. Equilibrium 

A* rt A rt 

response functions , T) , T may then be constructed through- 


(-> o 


.i* * * A* * def A * „* * 

(♦ » T1 » T ) - m (F , 6 ) = m(F , 0 , a(F , Q )) 


(7. 16) 


A* 


. * * Art 


A A A 


where m stands for any of f , T] , or T and m for any of ij/ 5 T] } .or T . 
Considering that 

• /\^ A /\ 

v 9 > = a e*& 0 > n > + e > n ) 3 e a <£> 0 > 


V~(F, 0) = 3_W O, n) + a J(F, 0, n)3 a(F, 0) 

r tl rv £ »*w 


It follows through the internal equilibrium equation that 


Art? rt * A * rt 

6 e** (£ ■ 6 > " V ( £ > 0 > ■> > 

«. A* . * _* . ^ * .* * 

V (£ . 6 ) ■ 5 f»(£ .«.»)■ 


(7.17) 


Therefore In. view of (6.11) and (6.12) 


* . A* * .* 

i) = - V* C£ . 0 ) 

* Art * 

& * V <£ • 0 > 

/■W 


(7.18) 


whe re 


* 1 * * T -1 
S “ ^ T ( P ) « 

#N° p* ~0 N /vr y 


(7.19) 


Equation (7.18) defines the equilibrium entropy and stress relations. 
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8. MATERIAL FRAME- INDIFFERENCE 

The constitutive equations are further restricted by the. axiom of 
material-frame indifference. This axiom states that the constitutive 
response functions must be form invariant under a change of frame. The 
change of . frame is characterized by a time-dependent orthogonal tensor Q(t) 
(cf. Truesdell and Noll [ 9 ]). 

Under a change of frame scalars 0, £, T| } and are unaffected. However 
Fj F, E, q j and T transform as follows 

F "-Q F 

*N> rsJ 

F - Q F 4- Q F 

fy rv rw 

g " Q 8 (8.1) 

fW <v 


9 "* Q 9 • 

rls 


T 




Q T 

r<i a^S 



The manner in which the porosity n transforms needs to be specified. It is 
postulated that n is unaffected under a change of frame. Therefore 

n ~ n, n -> n (8.2) 


under a change of frame. 

In view of .(,3.7) dependence on F can be indicated through dependence 
on L. Moreover L can be wrj tten as the sum of its symmetric part D and 

/v 

skew symmetric part W, 


L = D + W , 


T T 

where 2D = L + L , 2W *= L — L . Under a change of frame 

rv ^ rsv ^v> 


D 


QDQ 

/V fV» 


W 


Q W Q 

a; (V 


■ T 

+ Q Q . 


(8.3) 


(8.4) 
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• t *T 

Note that from the orthogonality of Q, it follows that Q Q = -Q Q . The 

rw AW aw Aw 

fact that W does not transform as a tensor under a change of frame can be 

Aw 

used to show that the constitutive equations should be independent of W. 

aw 

Therefore 


(T , q , n) = a(F, D, 6, g, n) 

AwS AwS Aw A* r Aw 


(8.5) 


/v A A A 

where a stands for either T, q, or n. 

AW Aw 

Necessary and sufficient conditions that the constitutive equations 
(6.10) and (8.3) satisfy the axiom of material frame indifference are the 
following 

a(F, 0, n) ="a(£F, 0, n) 

S.% <1* £' e ' &• ">S T - £<2 £■ SS2 T - e - 26- "> 


2 A <£■ £• e - b »> ‘ a<2b SS2 T - e - 2^. "> 


n(£, £> e > S. ■>) * n(Q £. wj , 6, Qg, n) 


( 8 . 6 ) 


A A 

where a stands for either f or T|. Choosing Q = -1, (8.6)^ ^ becomes 


AW 


T(F,.D, 0 > £’ n > = T(-F, D, 9, -g, n) 

Aw />*, aw ^ aw Aw 


-jfe £> e > & n > = <<-£» 9 > n > 


(8.7) 


A 

n 


(F, p. 


g, n) « n(-F j D, 0, 

rw «s/ »v 


n) . 


Thus T, n are even functions of F and g and q is an odd function of F and g. 

rv Ay/ fv rv 

Using standard arguments (cf. Truesdell and Noll [9]) it can be shown 
that a set of reduced forms of the response functions which are frame- 
indifferent are 
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where 


a(iF } , 0, n) = a (£, 0, n) 


T(F, D, 0, g, n) = F T + (C, C, 0, F T g, n)F T 
q(F 5 D, 6, g, n) = F q + (C, C, 0, F T g, n) 

f**/ w t**/ #N t r*J r* f** r*t f>j 


n(F, D, 0, g, n) = n + (C, C, 0, F T g, n) 

IN* »“W rw *w — -J 


T 

C = FF 


Is the right Cauchy-Green tensor. To arrive at (8.8) the identity 

• T 

C = 2F DF has been used. 

Considering that 

sJ(F, 0, n) = 2F 3 V + (C, e 3 n) 


equation (6.12) can be written as 


T = 2pF3_f*'(C J 0, n)F T . 

^jO (j tsj 

rw 


Also 


iq = -3 0 r.(£ J 0, n) 

T = 2pF3f + (C, 0, n)F T + F T + (C, C, 0, F T g, n)F T 


q - F q + (C> C, 0, F T g, n) 

^i»s A/ riy» ^ /%; ri 


n = n + (C, C } 0, F T g, n) 

rv> rw <SJ /C; 


and the equilibrium entropy and stress relations (7*18) become 


* , * 
f\ = - a flJ (C , 


_ A 

8 ) 


-f 

* ^ * *_ ,Vf ^ * 

£> - 2P l \J <£ . 



( 8 . 8 ) 


(8.9) 


( 8 . 10 ) 


( 8 . 11 ) 


( 8 . 12 ) 
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9 . MATERIAL SYMMETRY 


The constitutive response functions for the soil mass continuum' are 

further restricted by material symmetries. Recall Noll's L 10 3 definition 

of the isotropy groups?/ of a material response function F: the local 

isotropy group -%f of a material is the set of mass density-preserving changes 

A 

of local reference configuration which leaves the response function F 
unaltered. 


The above definition was arrived at through the recognition that a 
change from a given reference configuration to another reference con- 
figuration X which is indistinguishable from x by relating the values F 

. rwl 

A ' 

of F to deformation must be obtained by a mapping from X to X such that 



(9.1) 


For the soil mass continuum it appears more natural to base the defini- 
tion of the isotropy group $ of the soil for a constitutive response F on 
the solid aggregate mass density p and porosity n. Thus (9.1) is replaced by 






(9.2) 


or in view of (2.5) and (2.6) by the equivalent statement 



c 



(9.3) 


Therefore the definition of the isotropy group for a response function. F of 
the soil mass continuum is essentially the same as Noll's. Thus the isotropy 
group-2/ for the soil mass continuum is the set of all unimodular tensors H 
such that the following identities hold; 


/ a(l\ 0 , n) = a(FR, 0, n) 

T(F, C, 0, g, n) « T(FK, H T CH, 0, g, n) 


sO- 6, e. b n) - 2<rji, b t ch, e, j|. „) 
h b "> - S, T Si> e - b ■*> 


(9.4) 



23 


where a stands for either f or t]. 

It can be reasoned that if a solid is thought of as a body which 
has some preferred configuration from which any change of shape will change 
some of its properties then a soil mass is a solid. The ideas set forth 
by Noll [lO] can then be carried over unchanged, lie., since the soil mass 
continuum is a solid its isotropy groups is a subgroup of the orthogonal 
group $’(£) of linear transformations of Euclidean 3-spacc onto itself, 

M c (9(0 (9.5) 

When the soil mass-continuum is isotropic the isotropy group $ is equal to 
the orthogonal group &(£) . It can then be shown that the constitutive 
equations for fi, n, T^, q g> and n for an isotropic soil mass which are 
frame indifferent are 

♦ = ^(B, S, n) 

■T 1 = rf(B, 0, n) = “ Vl + (B, 0 ’ - n) 

T s = 2pa B (r + (B, e, n)B + if <B, B, 8, £, n) (9.6) 

q^ = q^(B, b, 8, g, n)- 

r&S >v fs* & 

n = n (B, B, 9, g, n) 

_ T • 

where B - F F is the right Cauchy-Green deformation tensor and all the 

response functions are isotropic tensor functions. Representation theorems 

for the response functions are available in [ 9 ]. 
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10. ELASTIC BEHAVIOR 


In solid mechanics elastic deformations are recoverable deformations 
and the thermodynamic process under which these take place -is a reversible 
process. It follows that elastic behavior occurs under' thermodynamic . 
equilibrium. Therefore for a soil mass continuum the constitutive equations 
for elastic behavior are 


tjf » lji*(F, 0) 

T * pO*(F, 0)F T 

fwS F ** 

q *= 0 


( 10 . 1 ) 


and 

n = d(F, 0 ) (10.2) 

wherein all deformations are recoverable and the temperature field is 
homogeneous. Note that (10.2) can be written as 


X(F, 0 , n) » 0. 


(10.3) 


Upon reduction of (10.1), (10.3) for material frame-indifference the 
following is obtained 

i)f * f^C, 9) 

1) » 0) 

Is ’ 2P 1 V'te e >f (10 - 4) 

q - 0 

X(£, 0, n) =0. 


For an isotropic soil mass (10,4) reduces to 
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♦ = t + (B f 0) 

>W 

T1 = -5 0 ^ + (B, 6) 

T = 2 P3 =, 0)B (10.5) 

*wS B V* n ' 


X(jB, 6, n) = 0. 


Inspection of (10.4) ^ ^ indicates that the elastic behavior of a soil 
mass is described by constitutive equations for i|f, T], T^, which are of 
the same form as the equations for a hyperelastic solid under isothermal 
conditions. However this is only possible whenever (10. 4)^ is satisfied. 
Consequently elastic behavior of the soil mass implies the existence X of 
(0,-0, n) such that (10.5),- is satisfied. 


11. LINEARIZATION 

In this section the constitutive equations postulated in the previous 
sections will be linearized by considering small departures from thermo- 
dynamic equilibrium. To this end it is assumed that the soil mass in its 
reference configuration is in a state of thermodynamic equilibrium. 

Introduce the Green-St. Venant strain tensor E which is related to C 
through 


' 2E = C - 1. 

^ rw 


(- 11 . 1 ) 


Consequently 

h = 0, n) = + 2E, 0 , n). 


( 11 . 2 ) 


Thus (8.11), (8.12) can be written as 
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*1 = -^(e, 0, n) 


T 

<"vS 




<fc+ 


= p F a'f(E, e,-n)F A + F T + (E, E, 0, F X g, n)F' 

** Cj a*# ^ awJj Vw ^ ^ ^ 

- F ^ + (E, E, e, F T g, n) 

A* ~ >W A»> J i*w A* ' 


T 

T.’ ^ 


(11.3) 


* A4- 

n = n (E, E, 

ru 


0, F g, n) 


• • 

where 2E = C has been used. 

Let be t ^ 1G displacement vector field of the soil mass. - The displace- 
ment gradient is then given by 


H = 


Vu(X, t) . 


rv r^J 


(11.4) 


H is related to the deformation gradient through 

F = 1 + H 

^ A*l 


(11.5) 


The departures from thermodynamic equilibrium can be measured by the 
quantity 6 defined by 


6 2 a (0-0 ) 2 + (n-n ) 2 + g • g 4- tr H H T + tr 11 H 1 , 

° O Ki K) <«w rv rj 1 


(11.6) 


where 0^ is the temperature in the reference state. The departure from thermo- 
dynamic equilibrium is said to be small if 6 < 1. A quantity of order S^is 
any scalar, vector, or tensor', denoted by 0( d ' /75 j, with the property that 
there exists a real number N such that 

II 0( 5 "^ I I = (11.7) 

as 6 -* o. Observe that 0(6 l/ ^)0(6 J7! -i) = 0 ( + . 

Under the assumption of 6 < 1 it follows that 
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0-0 — 0 ( 6 ) 
o 

n - n Q = 0(6) 

. H = 0(6) 
H = 0(6) 
£ = 0(6). 


Also 

E = ^ 4- 0(6 2 ) 

E = If 4- 0(6 2 ) 

F T g * g + 0(6 2 ) 

Av> rw 

|det £1 = 1 + tr |f + 0(6 2 ) 

p = P q (1 - tr |)+ 0(6 2 ) = P q 4- 0(6), 

where £ is the linear, strain tensor given by (3.14). 
Consider the ordered sextuple 

<£> i e > ■» 


( 11 . 8 ) 


(11.9) 


Its value in the reference configuration is 

(0, 0, 0 , 0, n ). 

V- f^/ o <%< o 

A+ 

Therefore f can be expanded about E=0, 8=0, n=n to yield an 

<v o O 

expression of the form 

1^(0 -!■ E , 0 4-0, n + a) = Ijf + a "0 4- \a^ 2 4- a.a'O + %a . a 2 + b, K] 9 

~ ~ o o oj. 2 3 4 1~ 


4- b_ C^T3ct 4 tr £E 4- k b-[l!f, El 4- 0(6 3 ) 


( 11 . 10 ) 
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where 


0=9-0 , a=n-n 

o o 


(11.11) 


and i{i , a , a , a are constants, £ is a second order tensor, b C * 3 > -b_[*] 

are linear and b„[* '] bilinear functions of L The constant tjf is the value 
3 ~ o 

of the free energy in the reference configuration. A term of the form a,_cc 

would ordinarily appear in (11.10) but because of the internal equilibrium 

3 

equation (7.14), = 0. Thus to within 0(6 ) 


f o t + aSB + IjaJ) 2 + a ,00 + %a.a 2 + b. [If] "0 + b^E] a + tr £ ? + %b,[Ef,E]. 
o 1 2 3 4 1 ~ 2 ~ ~ ~ 


( 11 . 12 ) 


It follows from (11.12) that approximate expressions for T|, for small 
departures from equilibrium are 


T| = -a.. - a.!? - a„a - b.. Ce] 
12 3 1 «— 


T 



m s + P 

~ ° 5e 


[S'] a + P £ + %p 

/•w 


O 






(11.13) 


Also the linear approximations- for small departures, from equilibrium of 

V " are 


T_. =' (b-TJ + b a)l + M. {S2 + M-Ctf] 
M) 1 Z ** '■vi 


q = "0 a + a a + K. [S'] + K.Cg] 

r^jS *vJL 'VZ ^ 


n = c-U + c a + N-Cfi*] + N 0 C8G 
1 2 1 ~ 2 ~ 


(11.14) 


where b. , b , c- c are constants, a., a„ are constant vectors, M, [•], 

± /. L Z ''-'1 *v2 aw 1 3 

M,[ - ] are symmetric linear tensor functions, K, [•], K_[*] are linear vccLor 

functions, and N^[- ], t * ] are linear scalar functions. In writing (11.4) 
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use was made of the fact that T_, n are even functions of g and q is an 
, rwD '"'•'S 

even function of E. 

The case of an isotropic soil mass is now considered- In this case 
b_ [• ] j b [•], ■£, b [• , ■ ] in (11.12) have the following representation: 

A A -J 


b LIT] = d tr E, a = 1,2 
■ a, ~ a ~ 

£ - d, 1 (11.15) 

b,K,ri - d.cr E 2 + d,(tr |) 2 

Also since there are no isotropic tensors of odd rank, a, , b ^, and [•] 
must drop out of (11. 4)^* and K^’]) M^[* ], N^C*], have the 

following representation: 

£2^ = K £ 


M, Ce*] = e. (tr-E)l + e 0 
M 2 m = e 3 (tr E)1 + e 4 

N, fif] = c tr S' ' 

1 ^ ^ 


N o m = c tr j£\ 

2 ~ 4 ~ 



ir 


(11.16) 


In (11.15), (11.16), d^, d^, d^, d^, d^, K, e^, e^, e^, e^, c^, and c^ are 
material constants. 

With the representations (11.15), (11.16) the linear constitutive 
equations (11.12), (11.13), and (11.14) become 

<= f + a^ + ^a/5 2 + a„CtT5 d.U tr E + d a tr Ii + d tr 'S' H- ljd.tr E 2 
of 2 3 1 ~ 2 ~ J ~ r*j 

+ %d,(tr E) 2 (11.17) 

_) ^ 

[Equation (11.17) continued on next page] 
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ri ® -a^aj -a , Ct - d . tr 

I 4 J 1 ^ 

* • 

T *= pi + (f/H + f„a + Xtr E + e tr 1)1 + 2\i 1 + c, E 
r^s ~ 1 2 ~ 3 rs/fv ~ 4 <*'•' 


JSs ” K & 

* 

n = cj + c.a + c tr E* + c trf 
1 2 3 ~ 4 ~ 


where 


P “ PoS > h ' P A + b l ■ f 2 " P o d 2 + b 2 

X = p d, + e, , 2\x = p d + .e . (11.18) 

OH- X O J Z 


The. genera] dissipation inequality (6.13) will give some inequalities for 
the constants b^, b^j c^, and the constants appearing in (11.16). 
However since the constitutive equations examined here arc* merely for 
illustrative purposes and have no bearing on true soil behavior these 
inequalities will not be found. 


12. -CONCLUSIONS 

In this study the basic equations of continuum mechanics for non-polar 
continuum have been reassessed as to their applicability to model a dry 
soil mass as a contimuum. It was found that the physical properties of a 
body which are taken in the construction of a continuum model, namely 
that of occupying a region of space and having mass, allow for the intro- 
duction of two densities, the bulk mass density and the solid aggregate 
volume density. In this respect the soil mass continuum differs from the 
strictly solid continuum. The introduction of the solid aggregate volume 
density brings the porosity of ‘the soil mass into play as an additional field 
variable. It must be pointed out that in applied soil mechanics the solid 
aggregate of a soil is assumed incompressible for computational convenience . 


for in this case the volumetric strain is computed through a phase diagram as 
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Av = A e 

v 1 + e ’ 
o 


( 12 . 1 ) 


e being the void ratio. Now it can be shown that in the limit the left 
hand side of (12.1) becomes tr E. Thus 



A e 
1 + e 


( 12 - 2 ) 


which was found earlier (see Eq. (3.20)). Now it is realized that soil is 
extremely difficult to sample and that volumetric measurements are hard to 
make. Nevertheless the assumption leading to (12.2) is dictated by this 
physical handicap and it is not in general a property of the soil. 

The equations used in continuum mechanics to depict the balance of 
mass, linear momentum, and energy also hold for the soil mass continuum 
provided that the mass density appearing therein be interpreted as the 
bulk mass density. In addition, through the balance of mass equation in 
terms of bulk mass density, a balance of mass equation in terms of bulk 
mass density, a balance of mass equation in terms of solid aggregate mass 
density was found. 

The introduction of the porosity as a field variable introduces 
complications since in general there is no equation relating porosity to 
deformation. This indicates that a constitutive equation for the porosity 
is needed. As an example of how a constitutive equation for the porosity 
may be introduced, a set of constitutive equations for a special kind of 
soil mass was studied. It must be pointed out that the constitutive 
equations studied here may not describe true soil behavior under load. 
However, the work presented here does yield results which are very important, 
namely that under isothermal conditions, 

(1) The constitutive equation for the effective stress should be of 
the form 

T =^(F,n) (12.3) 

r^G r+i <-w 

where indicates a general functional relationship. 
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(2) A constitutive equation has to be postulated for the rate of 
change of porosity. 

(3) Elastic behavior of the soil implies the existence of a function 
X such that' 

X(F,n) - 0 (12.4) 

for all elastic deformations. 
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SECTION II 

A THEORY OF SOIL PLASTICITY 
1 . INTRODUCTION 

A valid solution to a problem of the mechanical response of a soil 
mass to applied load must satisfy the basic balance equations of continuum 
mechanics. However, in order to obtain valid solutions, the constitutive 
relation between stress and strain of the soil must be established: Without 

knowledge of the stress- strain relation or its equivalent a so-called solu- 
tion is merely a guess. 

Partly because of the difficulty in obtaining self-consistent constitu- 
tive relations, problems in soil .mechanics are treated in several and 
unrelated ways. For example, -when no failure of the. soil is involved stresses 
at points in a soil mass under a footing, or behind a retaining wall are 
computed using linear elasticity.- Problems of bearing capacity, stability 
of slopes, failure of retaining walls are now being considered in the realm 
of plasticity, while settlement and consolidation problems are treated as 
essentially viscoelastic. 

In this paper our primary concern is with the establishing of a self- 
consistent phenomenological theory for the mechanical behavior of granular 
media which shows stress-strain behavior similar to that of an elastic-work 
hardening plastic metal. 

Drucker and Prager [ 8 ] suggested that the Mohr-Coulomb failure criteria 
for soils could serve as a yield function with which one could' associate a 
flow rule and to treat the soil mass as a perfectly plastic material. Even 
though several important and interesting results may be obtained by considering 
the soil mass to a perfectly plastic body with a' Mohr-Coulomb yield surface, 
predictions of volume changes, under this idealization, were higher than those 
found by experiments. Based on observations made on pressure-volume change 
curves Drucker, Gibson, and Henkel [7] subsequently explained that soil 
could be -treated as a work hardening material. Henkel [9], however, concluded 
that much of the available experimental, information for soils lay outside 
the scope of a useful theory of plasticity. Nevertheless Roscoe and co-workers 
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at Cambridge have indicated that certain soils can be described remarkably 
well by a simple isotropic work-hardening idealization- The work at 
Cambridge is discussed in considerable detail by Roscoe and Burland £ 5 ]. 

In the work described above a soil mass differs from a strictly solid 
mass by the constitutive relations describing the mechanical 'behavior . 
Aguirre-Ramirez [l] questioned the applicability of the basic balance 
equations of continuum mechanics to model a dry soil mass without investi- 
gation. He found that the two basic properties that bodies are known to have 
and which are used in the construction to model a body as a continuum, namely 
that (1) they have mass and (2) occupy regions of space, lead to the intro- 
duction of two densities one of which is the soil bulk mass density and the. . 
other can be related to the porosity of the soil mass. The basic local 
balance equations of a continuum can then be used for the physical- 
mathematical description of a dry soil mass and the processes occuring in 
it provided the mass density appearing therein be interpreted as the bulk 
mass density. The porosity, however, appears as an additional field 
variable that plays the role of a hidden variable. In soil mechanics 
porosity changes are related to bulk volume changes by an assumption which 
will be discussed in the text. 

In this paper we establish a self-consistent phenomenological theory.' 
for -the mechanical behavior of granular media which shows stress-strain 
behavior of the type discussed by Drucker," Gibson, and Henkel £ 7 ], This is 
done- by extending the Ideas presented in Cl], The main results of this paper 
are presented in Sections 4 and 5 we have felt compelled to include in 
Section 3 those aspects of soil mass behavior that we have used as a guide 
in arriving at them. In Section 4 we present a theory of plasticity for 
soils. The theory is a phenomenological theory -in which statements are made 
directly put into mathematical form and studied as such. Using the theory 
developed in Section 4 we construct constitutive relations for soils in a 
triaxial compression condition. This is done in Section 5. We- find that the 
Cambridge triaxial compression theory developed by Roscoe and co-workers £5] 
comes out as a special case of the theory constructed in Section 5. This is 
very promising because the Cambridge triaxial compression theory of Roscoe 
and co-workers [5] has been found to give reasonable agreement with experi- 
mental results. 

We must remark that by considering the soil mass to he dry we have 
disregarded the influence of the pore pressure on its behavior. This corre- 
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sponds to a soil in what is called a drained condition. 

NOTATION 

In this paper direct tensor notation is used except in Section 5. 

For the most part vectors in the three-dimensional inner product vector 

space U and points in Euclidean 3-space e are indicated by bold faced 

Latin minuscules: x,...,u. Linear transformations from U into U- are 

indicated by boldfaced Latin majuscules T,...,N. Second order tensors 

and linear transformations are regarded as the same. If T 'is a linear 
T -1 

transformation, T indicates its transpose, T its inverse, tr T its trace, 
and det T its determinant. The gradient with respect to spatial coordinates 
is denoted by grad and the gradient with respect to material coordinates . 
by V. 
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2. PRELIMINARIES 

We consider a dry soil mass body B which occupies a region R in Euclidean 
3-space in a reference configuration and denote by X the position in R of the 
particle X^B. We further suppose B to occupy the region R^ at time t and 
denote by x the position in R of the particle X^B. The motion of B from 
R to R fc is given by 

* = X(X,t>. (2.1) 

Let n^ = n Q (X) denote the porosity of B in R. According to the ideas set - 
forth by Aguirre-Ramirez [l] in order to describe 'the deformed state of B at 
time t we have to set alongside (2.1), 

n = n(X,t) (2.2) 

where n is the porosity of B at time t,‘ The function n is such that 

n o • n o<£) “ b(X,t 0 ) (2.3) 

where t is the reference time, 
o 

The gradient of 

I-'S') 

is called the deformation gradient. F is a second order non-singular 
tensor with the property 

i det f| > 0. (2.5) 

We let u = u(X,t) be the displacement vector from R to R and H be its 

~ r+j ~ t ^ 

gradient. 


H = V u(X, t) . 

The deformation gradient F is related to H by 

F = 1 + H. 

r\/ ^ Av 


( 2 . 6 ) 


(2.7) 
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We denote by dV and dv the element of solid aggregate of the soil body 
s s 

in R and respectively. The mean solid aggregate dilation or expansion A 
is then defined by [l] 


A = 


dv 

s 

dV ’ 
s 


It can be shown that A Is given by [l] 


( 2 . 8 ) 


A » 


O-Bl. 

(1-V 


Idct f| 


( 2 . 9 ) 


and is an additional strain measure that is characteristic of the soil mass 
body. 

The quantity e = e(X,t) defined by. 


n det £ 

0-n o ) 


( 2 . 10 ) 


i.s the void ratio of the soil mass defined as "the ratio of the element of 
void volume at time t to the element of solid aggregate of the soil mass in 
the reference configuration". Note 'that 

A = e> (2.11) 

n 


In what follows the word "specific" shall mean per unit mass of solid 

O O' 

aggregate. Let p , p be the soil bulk mass density and p , p the solid 

s s 

aggregate mass density in R and R respectively.. -The differential equations 
governing the deformation and motion of the soil mass body are given by [l] 

(i) Balance of mass 

pjdct F | = p° or p A = p°. ( 2 -12) 

~ s s 


(ii) Balance of linear and moment of momentum 


div T + pb = px 
T 

T = T 


( 2 . 13 ) 


(2.14) 



(iii) Balance of energy 


pe = tr T L - div q + pr. (2.15) 

In (2 . 13) and (2.14), T is Terzaghi’s [2 ] effective stress which was defined 
in [l ] in such a way so as to obtain (2.13) from a global balance law, b, 
the specific body force density, jx the acceleration, £ the specific internal 
energy density, q the effective heat flux vector, r the specific heat source 
density, and L the velocity gradient which is related to Y by 

L = F F* 1 = grad x (x,t) (2.16) - 

A» V { 

Alongside- (2.12), (2.13), and (2.15) constitutive equations are needed for the 
soil-mass. The work presented in [l] indicates that these constitutive equations 
should be 

f - Y(F,6,n) 

(T,q) = *5(F,F,e,g,n) (2J 17) 

n = H(F,F,0,g,n) 

'■v r*> r-* 

where ¥,$,5 indicate a general functional relationship ^ is the specific free 
energy density, 8 the temperature and. 

g = grad 0 . (2.18) 

^ 

Also in t 1 ] a triplet (F ,8 ,n ) with 

«i(F*,O,0*,O,n*) = 0 ' (2.19) 

was called a thermodynamic equalibrium state for the material point X of the 
soil mass. 

If the stress is written as the sum of a non-dissipative part T q and a 
dissipativd part T^, 
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with constitutive equations 


2o “£■><*> 6 > n > 

J D = &<£>£■ 0 >& n ) 


( 2 . 21 ) 

(2.22) 


where X > £~ indicate a general functional relationship, then £ must be 

r^j o ~D 
such that 


Zn (F*,0,6*,0,n*) = 0. 


(2.23) 


Also 

q “ 3C(F*,0, 9*,0,n*) = 0, (2.24) 

pu ^ /V <V* ^ 

i.e., at equilibrium the effective heat flux vanishes. In (2.24), K 
indicates a general functional relationship. 

In continuum mechanics elastic deformations are recoverable deformations 
and the thermodynamic process under which these take place is a reversible 
process. Therefore one may reason that elastic behavior occurs under thermo- 
dunamic equilibrium. Suppose the soil mass is responding elastically with, 
respect to some configuration R q at time t. Then according to (2.19) all 
deformations F, all porosities n, and all temperatures 0 are such that 

E(F,e,n) - H(F,O,0,O,n) = 0 (2.25) 

and since g = 0, the temperature field 0 is homogeneous. The effective 

/v • 

stress, free energy, and heat flux is given by 

T « T(F, 9, n) 

l|f * ?(F,9,n) (2.26) 

q = 0 
-- 

a *7 

where T,y are ordinary functions. 


+ r' 

Here a configuration and the region the soil mass B occupies in c. in 
that configuration are taken to mean the same. No confusion need arise since 
B is isomorphic to regions in £. 
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To demonstrate that the above Ideas of elastic behavior of the soil mass 
are equivalent to the elastic behavior (constitutive equation wise) of a 
simple continuum it suffices to indicate that the constitutive equations 
for an elastic soil mass are (2.25) and (2.26). Then under sufficient 

A 

smoothness assumptions on H the following equation 

n = f (F, 9) (2.27) 

for the porosity may be obtained from (2.25). It follows that 

T = T + (F,0) = T(F,0,f(F,0)) 

~ ~ ~ (2.28) 

♦ “ * + (F»0) = ?(F,e,f(F,0)) 

<v ' — ' 

which are the constitutive equations for an elastic simple continuum. 
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3. SOIL BEHAVIOR 

In this section we shall discuss some aspects of soil behavior which we 
shall use as a guide in the next section. Before doing this, however, let 
us indicate that most experimental data on the response of a soil mass under, 
load is for infinitesimal deformations under isothermal conditions. Also in 
the reduction of soil test data the following formula is used to compute 
volumetric strains 


Av _ An 

v 1-n 

o 


(3.1) 


where Av,An are the change in sample volume and porosity respectively. It 
can be- shown that under infinitesimal deformations (3.1) follows from (2.9) 
under the assumption of incompres’sibilily of the solid aggregate of the soil. 
In compressibility of the solid- aggregate is really not a general property 
of the soil. The reason for using (3.1) is that volumetric measurements on 
a soil sample are extermely difficult to make. Measurements of porosity 
changes on the other hand are simpler to make since these can be related to 
the amount of fluid expelled from the pores during the test. 

Most soils show a phenomena that is generally not shown by metallic 
solids undergoing infinitesimal deformations. This is the phenomena of 
dilatancy, i.e., bulk volume changes in a state of apparent pure shear. 

. Current methods of testing soils have been, almost always, restricted to 
the conventional triaxial compression test, the oedometer and to a far less 
extent, the direct shear test. The reader is. referred to Lamb and Whitman [ 2 ] 
for a .discussion of these tests. In the direct shear test the sample is 
assumed to be subjected to plane strain. In the triaxial compression test 
and ocdoincter the sample is assumed to be in a stress state in which the 
intermediate and minor physical components of the stress tensor are equal. 

The triaxial test is essentially a cylindrical sample first put under 
an equal all around pressure, called the confining pressure (denoted by o 
in Fig. 1), and then adding increments of load in the direction of the axis 
of the cylinder (Fig 1). As mentioned above the stress state of the sample 
during the test is assumed to be such that 



(3.2) 
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are the only non- vanishing physical components of the stress tensor in a 
polar coordinate system. In soil mechanics it is a practice to use the 
generalized stress parameters p and q defined by 

P = ~(a 3 + 2a) , q - a 3 - a ' (3.3) 

as an appropriate set of independent parameters that can be used in analyzing 
data obtained in the triaxial test. 

A typical stress-strain curve for a dry soil metss obtained in a vacuum 
triaxial test is shown in Fig. 2. This curve was obtained at the Geotechnical 
Laboratory of NASA's Marshall Space Flight Center at Huntsville for a lunar 
soil simulant material. For the purpose of discussion of this curve and for 
the remainder of the paper compressive stresses will be taken as positive. 
Examination of the curve (Fig. 2) indicates that the stress-strain behavior 
of this particular soil mass is similar to the behavior of an elastic-strain 
hardening plastic material. We also note that the response to a decrease in 
stress is an elastic recovery. 

The confined compression (oedometer) test is a cylindrical sample sub- 
jected to axial load but prevented from horizontal movement. Because of this 
last constraint lateral stresses develop which in general are not measured. 

As mentioned above accurate measurements of volume' changes in dry soil are 
not easy to make. In the oedometer test, however, .because of the no lateral 
movement constraint, the axial strain is exactly equal to the volumetric 
strain. The parameters used to analyze data from this test are generally 
the vertical effective stress denoted by p and the porosity n. 

A typical p-n curve obtained in the confined compression test is shown 
in Fig. 3. The significant features of this curve are the non-linear 
relationship between p and n and the elastic response to a decrease in stress. 
We note that ‘this curve also shows behavior which is similar to the behavior 
of an elastic-strain hardening material. 

The similarity of soil stress-strain curves of the type shox«n in Figs. 

2 and 3 to that of an elastic-strain hardening plastic material has led some 

+ 

researchers to suggest that the Mohr-Coulomb criteria could serve as a 

In the next section wc shall discuss the Mohr-Coulomb failure criteria. 
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yield function with which one could properly associate a flow rule [7], 
While this is a valid assumption it is not a reasonable one, because if the 
Mohr'-Coulomb envelope is used as a yield surface then yielding does not 
occur until failure takes place. As shown by the curve (Fig. 2) soils 
yield long before they fail. The use of the Mohr-Coulomb failure criteria 
as a yield surface -gives' erroneous predictions of high rates of change of 
volume during shear distortion. This is very unfortunate because research 
workers who reject these predictions will have the tendency to discount the 
usefulness of the theory of plasticity to model some aspects of the stress- 
strain behavior of soils. 
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4. SOIL PLASTICITY 

In this section we shall formulate a theory of plasticity for soils. 
While the considerations of Section 2 hold for finite deformations, the 
discussion in this and the following sections is limited to small strains. 
Thus the strain tensor that we shall use will be the linear strain tensor, 

E = %(H + H T ) . (4.1) 


Also we shall assume that the change in porosity 

C = n-n o 


(4.2) 


is of the same order of magnitude as the strains. 

In a soil we may define yield as a permanent irrecoverable deformation. 
We may write the strain E as the sum of an elastic or recoverable part 
and a plastic or irrecoverable part E^j 

E = E + E . (4.3) 

p 


In- addition we shall assume porosity to be given as the sum of an elastic 
or recoverable part and a plastic or irrecoverable part n , 


n = n + n , 
r p 


.(4.4) 


In view of this assumed resolution we shall have for the change in porosity 


c = 


+ n 

P 


(4.5) 


where C r = n r “ n 0 * The quantity v given by 

V = A - 1 


(4.6) 


where- A is given by (2.9) is the solid aggregate dilatation. For small 
strains V is given by 


V = 


(4.7) 
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Under the assumed ■ resolution (4.3), (4.5), v is given by 


C n 

v e, (_J_ + tr E ) + (-q — 2— + tr E ) . 
v l-n ^r y N l-n 

o o 


(4.8) 


Thus the solid. aggregate dilatation is the sum of- two parts, a recoverable 
part 

c 


V = + tr E 

r 1-n ~r 

o 


(4.9) 


and an irrecoverable part 


n 


V = + tr E . 

P 1-n. ~P 


(4.10) 


The plastic deformation of the soil is therefore described by the pair 

(E ,n ) . 

~P P 

A fundamental assumption of metal plasticity is that the elastic strains 
may be computed through the elastic constitutive equations for the stress. 

We carry this assumption into soil plasticity. Thus the elastic strains and 
elastic porosity may be computed through the soil elastic constitutive equation 
for an isothermal process. 


T - T(E C ) 

r* ^ i 


* 0 


(4.11) 

(4.12) 


where we have introduced the change in recoverable porosity. We assume (4.11) 


to be invertible in E , 

~r 


E„ » E(T,CJ. 

r* L ^ *** L 


(4-13) 


This allows us to write (4. 12) in the form 

2(t,£J - o. 


A. L 


(4.14) 
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We find it convenient to decompose the stress tensor into the sum of 
the mean pressure p and deviatoric stress i.e.. 


T = p 1 + 'f 


(4.15) 


p = -r tr T, tr 


T= o. 


( 4 : 16 ) 


In view of this assumed resolution we may write (4.14) as 


g(p,T>C r ) - o 


(4.17) 


Equation (4.17) defines a surface fi in (-T space. We call this surface Q, 
the elastic surface. Note that we may solve (4.17) for so as to obtain 


g rf (p/n. 

A/ 


(4.18) 


The curve in C-p space defined by 

C r * = C r (p) = g*(p,0) (4.1-9) 

is called the elastic swelling curve of the consolidation curve of the soil. 
We consider a curve in £-p space given by 

♦fp,C) = 0 (4.20) 

with fixed and unique. For a given soil there exists such a curve which 
is called the virgin consolidation curve . The intersection of the elastic 
swelling curve with the virgin consolidation curve defines a point in C-p 
space which is a yield point for the soil. We denote the mean pressure 
corresponding to such a point by p Q . Now (4.20) may be solved for the change 
in porosity 

C = C(p) . (4.21) 


From (4.19) and (4.21) we find 


n p = 


£r ( P 0 > 


(4.22) 
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which indicates that with the pressure at yield we can associate a plastic 

porosity n . We shall use this result below. 

P 

One of the main ingredients of metal plasticity is that of a yield surface. 
We shall now demonstrate how to construct a yield surface for a soil. Our 
construction of a yield surface for the soil is wholly dependent on the 
hypothesis that "plastic porosity is a unique function of stress,” i.e., 

n = N(T) . ’ (4.23) 

P 

We further assume (4.23) to be the solution for n^ of the equation 

G(T,n ) = 0 (4.24) 

~ p 

where G is unique. In plastic pososity-stress space (4.24) defines a six- 
dimensional hypersurface S that is called' the state boundary surface. A 
soil particle will be said to be in a plastic state if the value of the 
stress and plastic porosity are such that (4.24) is- satisfied. 

We consider a curve on E. The projection of this curve on stress space 

is a five-dimensional hypersurface Y. There are curves on E which have the 

unique feature that the value of plastic porosity is the same all along the 
* 

curve. Let n^ be the fixed value of the plastic porosity along one of these 
curves and consider the set B of all T such that 

G(T,n*) = 0. " (4.25) 

r 

We call B a yield domain. The projection of this equi-plastic porosity 
curve on stress space is a five-dimensional hypersurface whose equation is 
given by 


A(T,x.) = 0 , T£B 


(4.26) 


where the parameter H depends upon the value n^ of plastic porosity, i.e., 

A * 

H = X(n ) . 

P 


(4.27) 
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Our concept of a yield surface Y, for the soil is given by (4.26). Note 

that along the- yield surface the plastic porosity has a constant value 

n , i.e. , for all states of stress which locate points on Y the value of 

the plastic porosity is the same. In particular we can find a yield 

surface such that the value of plastic porosity associated with it is given 
* o 

by (4.22), i.e., n^ = n^. We call such yield surfaces "volumetric yield 
surfaces" and denote them by ¥ . . Thus a volumetric yeild surface is 
characterized by 


A O . 

K = X(n ) . 
P 


(4.28) 


Also since with each n we can associate a mean pressure p we can also 

p ° 

characterize volumetric yield surfaces by 


H = K (p o ) 


(4.29) 


If we consider another equi-plastic porosity curve in- £ its projection 
in stress space is given by an equation of the form (4.26) but with a 
different yield domain and, of course, a different value of the parameter 
corresponding to a different value of plastic porosity. In addition, since 
for all states of stress on a given yield surface, the plastic porosity is 
constant, it follows that the change in porosity Q is also constant. We 
have envisioned soil as a work- hardening material. Thus in order for plastic 
deformation to occur, the stress point must move outside the yield curve, 
i.e., the initial yield point must be exceeded. A new yield curve is then 
established which, depending upon the shape of the state boundary surface £, 
may or may not resemble the old yield curve. We shall also show below that 
the yield surface for soils is not a closed surface. 

With our concept of a yield surface the loading and unloading criteria 
are respectively given by 

tr 9 m A(T,K)T < 0, A(T,k) = 0 

tr b A(T,h)T > 0, A(T,H) = 0 

1 <V 


(4.30) 
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while the neutral loading is given fay 

tr 3 A(T, ?QT = 0, A(T,K) « 0. (4.31) 

N 

An important' ingredient of plasticity is that of Drucker's postulate 
of stability of material [ 3 ]. This postulate is used to classify a material 
as a work-hardening material. According to this postulate, if an external 
agency applies a small surface fraction which alters the stress at each point 
by T, then upon gradual application and removal of this surface fraction 

tr T E * 0 (4.32) 

rsi «V>p 

if the material is work-hardening. Important consequences of Drucker's 
postulate are [ 3 ]: 

(i) The yield surface and all subsequent loading surfaces must be 
convex. 

(ii) The plastic strain increment vector must be normal to the loading 
surface at a regular point, and it must lie between adjancent normals to 
the loading surface at a corner of the surface. 

The normality condition (ii) implies that at smooth points on the 
yield surface 

E = X P(T,K) (4.33) 

is. p r~t ~ 

where \ is a function of the deformation history and is such that (4.33) is 
homogeneous in time and 

D(T,H) = 3 A(T,H). 

A- 1 ft- 

IV 

We can also write the yield surface in the form 

^(P,T^) = 0 


(4.34) 


(4.35) 


The normality condition then leads to 
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e = X D(p,°r>) 

^p f' 

V- " a (p/.^) 

P " 

where 

dCp.'Jk) = ^.ACp.r.n) 

A» 

a(P»°T» k ) = 9 A(p,T, K) 

" p 


(4.36) 


(4.37) 


and e is the plastic strain deviator, cp the plastic volumetric strain 

P 

given by 



1 

3 



cp «= tr E . tr e = 0. 

~p »** p * "p 


(4-38) 


We note that the dependence of A on'T must be such that 

tr D( P} T,k) = 0. 


(4.39) 


We shall now indicate how the Mohr-Coulomb failure criteria may be 
used together with our concept of the yield surface. The Mohr-Coulomb 
failure criteria states that the- magnitude' of the shearing stress' T on any 
section through a mass of an isotropic cohesive soil must not be greater than 
an amount which depends linearly upon the normal stress O acting on the 
section. This condition is expressed as 

T £ c + a tan $ (4.40) 

where c is the cohesive and $ the angle of friction of the soil. Failure 
can occur when the equality sign in (4.40) holds for some section through 
the soil. Shield [ 4 ] has constructed a surface in principal stress space 
corresponding to 


T = c + C7 tan $ 


(4.41) 



52 


Letting Cf cr^o^ be tbe principal stresses with CT^ 2 s s a^, shield E4 ] 
finds that the surface is a right hexagonal pyramid equally inclined to 
the a 3» a 2 ,C7 3 axes anc * w ^th vertex at the point = CT^ = < 2 ^ = -c cot § 

(Fig. 4). When the stress point of the soil is on this surface, the soil 
mass is said to be on a -failure condition. The significant feature of the 
failure surface is that it is not a closed surface in stress space but 
rather an open .surface. Consequently it divides stress space into two 
regions I and XX. Stress points in region I are such that the inequality 

< c + O' tan $• (4.42) 

is violated and consequently such states of stress are not possible for the 
soil since by definition the soil has failed. On the other hand, points in 
the region II are such that the inequality (4.42) holds and consequently 
equilibrium of the soil mass is possible. However, since soils yield 
before they fail, all points in region II in the vicinity of the failure 
surface also lie in some yield surface. For a soil which has work-hardened 
isotropically two possibilities arise: (1) ¥ is tangent to the failure 
surface, or (2) ¥ traverses the failure surface. If the possibility. (1) 
prevails the plastic increment vector at the tangent point will be normal 
to both surfaces. Considering that the use of normality of the plastic 
strain increment vector to the failure surface gives erroneous predictions 
of high rates of change of volume during shear distortion we disregard 
the possibility (1) above. Therefore the failure surface is traversed by 
the yield surface. This indicates that stress points can only lie on a por- 
tion of the yield surface, that portion which lies in region II. 

Consequently, the yield surface is not a closed surface. 

We now assume the soil to be isotropic. We further assume to be 
linear in T. Therefore (4.13) has the representation 

S r = < a i(£ r > tr P 1 + a 2 <C r )T (4.43) 

1 > 

where a^ and are functions of the change in recoverable porosity. Letting 

J_ = trf 2 , Jo « tr-T 3 , 


(4.44) 



(4.18) reduces to 

C r = g*(p, J 2 , J 3 ) (4.45) 

Also we assume the soil to be isotropic work-hardening material and use the 
form (4.35) of the yield surface. Thus 

ip.'T, H) = A(p, J 2 , J 3 , *) (4.46) 


From this follows 


D(p,7>) = A(p 3 J 2 ,J 3 ,K)r+ ,J 2 ,J 3 ,Kff 


(4.47) 


and cpnsequently 

tr = 3S j a(p,j 2 ,j 3 ,h)j 2 . 

In view of (4.39) and the fact that J„ ^ 0, this last equation leads to 


a j ^A(p,J 2) J 3 ,H) » 0 


(4.48) 


or equivalently that A must be independent of J^. Consequently (4.35) reduces 
to 


A 


A(P, J 2 , H) = 0. 


(4.49) 


Considering that 


j 2 - h - 3 p 2 • h ■ tr i 


(4.50) 


we can write (4.49) in the form 


A(p,J 2 ,K) = A(p,I 2 ,X) = 0 


(4.51) 


or using principal components of stress 
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We indicated at the beginning of the last section that (3.1) follows 
under the assumption of incompressibility of the solid aggregate. A formula 
similar to (3.1) may also be derived by recognizing that fundamentally there 
are two mechanisms that contribute to the deformation of the soil; 
distortion of individual particles and relative motion between particles as 
the result of sliding or rolling. If we assume that plastic deformations 
are mainly due to the relative motion of the soil particles due to sliding 
and rolling then a reasonable assumption we may make is that during this 
motion the soil particles are essentially incompressible. Mathematically 
this assumption is stated as 

V » 0 (4.53) 

P 


where V is the irrecoverable part of the solid aggregate dilatation given 
P 

by (4.10). In view of the ass mption (4.53) we obtain from (4.10) . 


n 



(4.54) 


which is a formula similar to (3.1). 

Now considering that 

A * A • . .. A - A . 

A = ap + tr D + 3 A 3 ftn = 0 
•v I ft n p 

p r 

A 

where o. and D are given by (4.37), we can write (4.54) in the form 

• . A • A 

e = A(ap + tr D T)D 

P +* A/ ^ 


cp = A(ap + tr D T )a 

p 


(4.55) 


where 


(1-n ) SJ 3 ft a ‘ 
v o' ft n 

P 


(4.56) 


Thus under the assumption (4.53) the flow rule is given by (4.55). 
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5. STRESS- STRAIN THEORY FOR TRIAXIAL COMPRESSION CONDITION 

In this section we shall construct from the theory presented in the 
previous section stress- s train relations for the independent interpretation 
of triaxial compression condition. With reference to Fig. 4 we assume, that 
in a triaxial compression condition, the non-vanishing physical components of 
the stress tensor to be 


a 


1 


J 



We further assume a homogeneous stress field in which case it follows .from 
equilibrium that 


O 


1 



(5.1) 


Under these assumptions the mean pressure and physical components of the 
deviatoric stress tensor q^ (i — 1,2,3) are 

P = i(c 3 + 2a 2 ) , , 3 - |(0 3 - op 


(5.2) 



If we assume an isotropic material then the principal directions of the 
stress tensor and strain tensor coincide. The non-vanishing physical com- 
ponents of the strain tensor are then given by 


£ 


1 




and the physical components of the strain deviator by 


e 




3 



e 


1 


■\e 


3 


(5.3) 


where cp i s the volumetric strain. 
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We can write the stress power of the soil mass in the form. 


D = tr T E = trT e + pep. 


(5.4) 


Using the principal components of and e we bring (5.4) into the form 


D = qe + pep 


(5.5) 


where 



(5.6) 


The form (5.5) suggests that in a triaxial compression condition the 
generalized stress and strain parameters be (q,p), (e,cp). Accordingly we 
decompose (e,cp) into the sum of recoverable and irrecoverable parts. 


e = e + e 
r P 

(5.7) 

cp = cp r + cp ? 


and construct a yield function in q-p space. 


In two dimensional stress space the yield surface reduces to a curve 
and the Mohr-Coulomb failure surface to two straight lines meeting at a point 
on the p-axis. In Fig. 5 we show these two lines for a cohesionless soil 
(c “ 0) . The angles ^ are defined by 


tan 


6 sin i 
3 - sin § ’ 


tan $ 2 


6 sin $ 

3 + sin $ 


(5.8) 


where $ is the angle of internal friction of the soil. We note that the two 
Coulomb lines 0C^, OC^ divide q-p space into two regions I and II. In Fig. 5 
the current yield locus Y is the curve F^F^ w ^- c ^ traverses the p-axis at p Q . 
Recall that along the yield curve the value of the change in porosity is 
constant. This value of the change in porosity can be obtained from the 
soil virgin consolidation line which is obtained through another experiment. 
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Without loss of generality we take 
less quantities 



and introduce the dimension- 


M = tan 

N c (5.9) 

P o 

11 a ^ ( 

P 

The value of p,q at the point is then given by 

q = Mp (5. 10) 


where' a bar over a quantity indicates that these are the values of q,p when 
the soil is in a failure condition. We also introduce the dimensionless 
quantity N£ defined by 


N 


f 


2 

Po 


(5.11) 


We note that whereas M is a constant for the soil, is constant only for 
a particular yield. However for lack of experimental evidence to use as a 
guide we shall consider N, to be a constant for the soil. This assumption 
in itself suggests experimentation. 

We shall assume the current and subsequent yield locii to be symmetrical 
about the p-axis and to be segments of ellipses which pass through the origin 
It can then be shown that the equation of the ellipse is 

N= « ( 5. 12) 

M + KX| 


where 




(5.13) 
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Thus the current yield locus is given by (5.12). In (5.12), p Q appears as. 

a parameter that we may use as a strain-hardening parameter. This being the 

case the yield locus will be a volumetric yield locus. The dependence of 

p on n has to be determined in order to obtain a flow rule associated with 
o • p 

(5.12). To this purpose we assume that the virgin consolidation line of, the 
soil when plotted on Q - In p space is given by Terzaghi's well known equation 


P Q 

C « -X ln(— ) 

0 P ± 


(5.14) 


where X Is a soil constant and p Is the consolidation pressure of the soil 
o I 

in its reference state. We also assume the elastic swelling curves to be 
straight and parallel lines of slope k when plotted on £ln P space. We can 
then show that the construction depicted by (4.22) leads to the following 
relation 


-n 

P q = Y exp(~~y£) 


(5.15) 


where Y is the initial yield pressure under confined consolidation and 


Y ~ X - K . 
o 


(5.16) 


In order to construct a flow rule associated with (5-12) we assume that 
(4.54) holds. We can then show that the coefficient of proportionality X 
in (4.36) is given by 


x- = 


Yp 

— + 
(1-0 p 


2KT|n ) ( M i 2 + KTj 2 ) 


2 2 
M + Kt| 


M 


KT) 


(5.17) 


The plastic strain rates can then be shown to be given by 


; =— Y fJELwP+JSL) 

P M 2 - KT ) 2 p M 2 + KT ) 2 


i = Y ,P 2KT1T1 

( P + kt, 2 ’ ' 


(5.18) 
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Having obtained expressions for the plastic strain rates the next order of 
business is to obtain expressions for the recoverable strain rates. It can 
easily be shown Lhat (4.43) and (4.45) reduce to 

e = )q 

r 3 r ^ 

<p r = 3 V C r )p (5 ' 19 > 

G r = 8 (P,q) 

where 

a 3 (C r } ’ V C r> " 3 W + a 2 (C r ) (5 ‘ 20) 

For lack of experimental information regarding the dependence of C r on 
q we assume 

3 q g*(p s q) = 0 , (5 • 21) 

* 

i.e., g is independent of q. It follows that 

C r =g( P )- (5.22) 


This being the case, we have already assumed the form of (5.22). Equation 
(5.27) describes the elastic swelling curve. Consequently 



(5.23) 


Using (5.19)2. (5.22), and (5.23) we find 

¥ r = 


(5.24) 


where . 


t(C r ,p) - 3p(a 4 (C r ) - 3, a (C r )). 


(5.25) 
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If we assume that there is no recoverable energy during shear distortion 
then 


e = 0. 
r 


(5.26) 


The final stress-strain relations for triaxial compression condition 
are obtained by combining (5.18), (5.24), (5.26), and (5.7), 


where 


4 = 4 = V . f 2 KT 1 £ 2 KTTH , 

r (1-n ) ^ 2 2 ^d + 2 2^ 

' o' M - K Tf P M + K if 


i i r-gYK m , , + X * (C p) Pv 

1 P ( l -o> M 2 + Kt, 2 (?I ’ P) P ) 


(5.27) 


(C r ,p) = Y + «(l-n o )f(C r ,p) . 


(5.28) 


Let us now consider the following form of the material function a. (C ) 

4 r 

appearing in (5.19)2, 


V C r>. = 3R(l-n )(!-«) < exp C r ’ exp <-iT» 


(5.29) 


where R is the residual pressure, i.e,, the pressure experience by the soil 
when it has been held at rest in its reference configuration at all times. 
The assumed equation of the elastic swelling curve is 

-c 


p = R exp(— ) . 


(5.30) 


Substitution of (5.29) and (5.30) into (5.19) yields 


Cp = 

V r (l-n o )(l-K) 


“C ( 1 -K) 

(exp( — 1). 


(5.31) 


If we expand the exponential function into a power series we can write (5.31) 
in the form 

*r - wh + 0 O • 

0 
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But since we have assumed cp^_ and Q - to be of the same order of -magnitude, 
it. follows that 


" (l-n ) 


(5.32) 


which is an equation of the form (3.1). Therefore if we assume (5.29) we 
may replace (5.19)^ by (5.30) which may be written in terms of Cp^ as 

ln <R>’ ’ (5 ' 33) 

^ " * 

Under the assumed form (5.29) for a (£ ), the function X (G ,p) given by 

f ^ 

(5.28) reduces to 

X*(C r ,p) = X (5.34) 

r o 


In this case equation (5.27)2 reduces to 


cp = 


2YWI 


^""o^ M 2 + KT1 2 


+ x £) 


o p' 


(■5.35) 


while (5.27)^ remains the same. 

Let us assume the particular value of \ for the .coefficient given 
by (5.11). In this case K c 1 and the yield locus and stress-strain relations 
(5.27)^ and (5.35) reduce to 


N = 


M 


2 2 
M + T] 


i „-_X— ( — 2T] E . 2TW , 

( l-n ) \ 2 2 M p t 2 , „2 J 

v o^'M-Tl M + T) 


(5.36) 


cp 


+ x £> 

( l-n ) '2 L 2 + o p ; • 
v o' M + T| ^ 


# “j- 

If we assume the line 0C^ in Fig. 5 to be critical state line instead of a 


"*~See Schoffield and Wroth C6 ] for a thorough treatment of the critical 
state concept in soil mechanics. 



Coulomb line then (5.36) is identical to the Cambridge triaxial compression 
theory presented by Roscoe and Burland [5] for "wet" clay. We must point 
out that Roscoe and co-workers at Cambridge have, for the past decade, 
concentrated considerable effort to arrive at self-consistent constitutive 
relations for soils. The reader is referred to the recent article by 
Roscoe and Burland [5] for an account of the work at Cambridge. 
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SECTION III 

ON A YIELD SURFACE FOR SOILS 
1. INTRODUCTION 

In 1952, Drucker and Prager [ 3 ] introduced- an idealization for the. 
phenomenological behavior .of soils under load.. In this idealization the 
soil mass, is treated as., a perfectly plas.,. tic material 1 with the. Mohr- Coulomb 
failure .criteria for soils. as a yield function with-which a flow rule can 
be associated. Volume changes predicted under this idealization, however, 
were higher than those found by experiments. Based on observations made 
on pressure-volume change curves Drucker, Gibson, and Henkel [ 4 ] explained 
how soils could be treated as a work-hardening material. Roscoe and 
co-workers at Cambridge have indicated that certain soils can be described 
remarkably well by a simple isotropic work- hardening idealization: The 

work at Cambridge is discussed in considerable detail by Roscoe and 
Burland [ 2 ]. 

'Aguirre-Ramirez and Costes [l ] presented a self-consistent phenomeno- 
logical theory for the mechanical behavior of granular media which shows 
stress-strain behavior of the type discussed by Drucker, Gibson, and 
Henkel [ 4 ]. As an example a yield surface and associated flow rule, for 
triaxial compression conditions, was constructed in [l]. It was also 
shown in [l] that the Cambridge triaxial compression theory developed by 
Roscoe and co-workers [ 2 ] comes out as a special case of this theory. 

In this paper we generalize, to complex stress fields, the special 
triaxial compression theory presented in [l]. We have been encouraged to 
do so in view of the fact that the Cambridge triaxial compression theory 
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of Roscoe and co-workers has been found to give reasonable agreement with 
experimental results. Accordingly in Section 2 we briefly review the 
theory of soil plasticity presented in [l]. In Section 3 the special 
triaxial compression theory of [l] is also reviewed. This theory is gen- 
eralized to three-dimensional complex stress field in Section 4. Recog- 
nizing that many soil mechanic problems can be idealized to plane-strain 
situations we present in Section 5 a theory for plane- strain. In Section 6 
we compare both the three-dimensional and plane-strain theories to those 
of Roscoe and Bur land [ 2 ] and find perfect agreement. 
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2 . PRELIMINARIES 


.We consider the particles of a dry soil mass continuum to 
be referred to a fixed rectangular Cartesian system x_^ (i=l,2,3) 
and let u denote the displacement field. The. linear strain tensor 
is then given by 


Uj = V 2 t u i,3 


u j,i> 


( 2 . 1 ) 


where we have denoted partial differentiation with respect to x^ 
by ( ),£. We also let n,n Q denote the current and initial poro- 
sity of the soil mass continuum respectively and define the 
change in porosity £ by 


' £ = n - n Q . (2.2) 

We further assume £ to be of the same order of magnitude as the 
linear strain tensor y . . . 

Earlier we have defined yield in a soil to be a permanent 

irrecoverable deformation [1] .- The strain may then be written 

as the sum of an elastic or recoverable part y|j and a plastic 

or irrecoverable part y('., 

^•3 


13 


= Y 


ID 


+ Y 


ID 


(2.3) 


In addition we shall assume porosity to be given as. the sum of 
an elastic or recoverable part n’ and a plastic or irrecoverable 
part n", 


n = n ' + n" 


(2.4) 
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In view of this assumed resolution we shall have for the change 
in porosity 

£ = £' + n" (2.5) 


where £ ' = n ’ - n o . 

According to the ideas set forth in [1] the irrecoverable 
deformation of the soil is described by the pair (y^ 11 , n”) . 

The theory of soil plasticity presented in [1] is wholly dependent 
on the hypothesis that there exists a unique function F of plastic 
porosity n" and Terzaghi’s effective stress such that 

F(a- n") - 0. (2.6) 

X J 

In plastic porosity-stress space (2.6) defines a six-dimensional 
hypersurface i which was called in [1] , the state boundary sur- 
face. A soil particle was then said to be in a plastic state if 
the value of the stress and plastic porosity are such that (2.6) 
is satisfied. 

A yield surface ft for the soil was defined in [1] to be the 
projection on stress space of curves on £ along which the plastic 
porosity has the constant value n " . This is given by 


where 



k) 


= 0 


(2.7) 


k = k(n") (2.8) 

is the strain-hardening parameter. The flow rule associated 
with (2.7) is given by 
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y" = X d ( a , k) 
mr rnr ij 


where X is a function of the deformation history 
(2.9) is homogeneous in time. Also in (2.9) 


D (a . . , k) = 

^ ^ 3 a 


3A 


(o ^ j * k) . 


mr 


We often find it convenient to decompose 
of its deviatoric part and mean pressure p. 


cr . . = p6 • . + x • • 
a-3 p 3-D 3-D 


where 6, . is the Kronecker ' s delta and 
3-D 


p = 1/3 a t 
mm 


x = 0 
mm 


Under the resolution (2.11) of the sttess we may 
the form 


A.(p* x . . / k) — 0. 


The flow rule (2.9) then reads 


e" = X D (p/ x . . , k) 
mr mr ij 


Y" = X a(p, ± ., k) 
mr 3 


where e" is the plastic strain deviator and 
mr 

D mr (p ' T ij' k) “ T i y k> 


mr 


3A , 


a(p, x.., k) = — (p, x ii , k) 
13 3p 3 


. (2.9) 

and is such that 

( 2 . 10 ) 

, into the sum 

( 2 . 11 ) 

( 2 . 12 ) 

write (2.7) in 

(2.13) 

(2.14) 

(2.15) 


Fundamentally there are two mechanisms that contribute to 
the deformation of the soil: distortion of individual particles 
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and relative motion between particles as the result of sliding 
or rolling. We assume that soil irrecoverable deformations 
are mainly due to the relative motion of the soil particles due 
to sliding and rolling and that during this motion the soil 
particles are essentially incompressible. It can then be shown 
that under this assumption [ 1 ] 


(" = - 
'mm 


n' 


u-V 


( 2 . 16 ) 


The flow rule (2.14) may now be written as 


e" = A(D t . . + ap) D 
mr 13 31 mr 


(2.17) 


Y n 85 A(D t. + ap)a 
ittn 1 ! 


where 


. r / . v BA Bk . “i 
A - ■[ (l-n 0 ) a — ] 

Bk Bn” 


( 2 . 18 ) 


The recoverable strains y* and porosity n' may be computed 

mr 

through the soil elastic constitutive equations [ 1 ] 


Y 1 = Y* (o. , , V ) 
mr mr 1 j 


S' = V (o ij ) . 


(2.19) 


f c> . . 

For an isotorpic soil for which the strains are linear m a, , , 

^ J 

(2.19) j has the representation 


Y * = £ + a 2 U')a (2.20) 

mr 1 ii mr mr 

where ai, a 2 are material functions of recoverable porosity 
We may also write (2.20) in the form 
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e mr - (2 - 21) 

'C = 3 “» (t ’>P 

where 


a 3 U ') = 3ai U ' ) + a 2 (? *) . 


( 2 . 22 ) 


The following form of (2.19) was assumed in [1] 

£ 1 — ~K ln(£) (2.23) 

K 

where k is a soil constant and R is the residual pressure, i.e., 
the pressure experienced by the soil when it has been held at 
rest in its reference configuration at all times. We note that 
in view of (2.23) we may write (2.21) 2 in the form 


Y mm = 3Ra 3 U ') exp ( -— -) 


(2.24) 


which relates volumetric strains to recoverable change in poro- 
sity. In [1] it was shown that under a suitable chosen function 
a 3 the relation (2.24) can be reduced to 


y • = - 

'mm 


JLl 


(l-n 0 ) 


(2.25) 
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3 . TRI AXIAL THEORY 


In this section and for the duration of the paper we shall 

consider the soil mass to be isotropic. We shall also use the 

convention that compressive stresses will be taken as positive. 

We consider the family of deformations y such that when 

1 mr 

referred to a suitable set of orthogonal axes 


v = °' m * r 

Yu = y 22 - 


(3.1) 


Such deformations are called triaxial deformations. Since we 
have assumed the soil to be isotropic it follows that the prin- 
cipal directions of the’ strain and stress tensor coincide. This 
being the case then it follows that under triaxial deformations 
the components of the stress tensor take the particular form 


mr 


= 0 , m 7 ! r 


1 1 


2 2 


(3.2) 


It can be shown [1] that under triaxial deformations a suitable 
set of generalized stress and strain parameters are (q,p), is, 8 ) 
where 


*3 ~ ^33 ~ cr , . t r> “ ( O' 3 3 + 2cr j j ) 

* (3.3) 

e = 2/3 (y 3 3 - Yll ) , 0 *= Y^ = Y 33 + 2y a 1 . 

In terms of these generalized parameters, the stress power is 
given by 



77 


D = q e + p0 . (3.4) 

We consider the following invariant. of the deviatoric 
stress tensor 

J = t x (3.5) 

mr rm 

It can be shown that in a triaxial stress field j is given by 

J = 2/3 (a 3 3 - cxi i ) 2 = 2/3 q 2 . (3.6) 

The octahedral shear stress T 0 is related to-J through 

x 0 = . (3.7) 

3 

In a triaxial stress field (3.7) reduces to 

To = /2/9 q. (3.8) 

Equations (3.6), (3.8) give meaning to the stress parameter q 
in terms of invariants of the stress tensor. 

The strain measures e , 0 may be decomposed into the sum of 
recoverable and -irrecoverable parts 

e =' e ' + e" , 0 = 0 ' + 8 ". ’ (3.9) 

In .[1] a one-parameter family of yield curves was constructed 
in q-p space for a cohesionless soil with an angle of internal 
friction $ at failure. A member of this family is shown in Fig. 1 
and is given by 

_E = 

p 0 m 2 + Kn 2 


( 3 . 10 ) 



78 


where n = g/p and 

6 sin* 

M =» 

3 - sin<j) 

Here is the value of the mean pressure at failure at a change 
in porosity 5 and p 0 represents the pressure corresponding to 5 
on the virgin isotropic compression curve of the soil. The para- 
meter of the family was taken in [ 1 ] to be p 0 and is given by 

po = Y exp (-^-) (3.12) 

3 

where Y is the initial yield pressure under confined consolidation- 
and 

3 - Xo ~ K (3.13; 

A 0 being a soil constant. 

Under an assumed constant value of K the following flow 
rule was established in [ 1 ] , 

e" — § ( — 2Kt I .) (E + ^Kriri — j (3.14) 

(1-no) M 2 -Kri 2 p M 2 +Kr; 2 

■ q ' — _3. ( e + mi -) ; 

( 1 -no) p M 2 +Kn 2 

In [1] we indicated that the yield curve given by (3.10) • 

and associated flow rule (3.14) reduce to the "Cambridge" triaxial 
compression theory presented by Roscoe and Burland [2] if we 
interpret the line OC in Pig. 1 as a critical state line instead 
of a Coulomb line and if we take K - 1. 


K = 


Po - P u 


u 


(3.11) 
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. 4. THREE DIMENSIONAL THEORY 

The yield curve given by (3.10} and associated flow rule 
(3.14) is in terms of triaxial compression stress parameters 
q-p and as it stands it is only good for analysis of triaxial 
deformations under a triaxial stress field. In general, however, 
it is desirable to obtain a yield surface and associated flow 
rule to analyze deformations under complex stress fields. Now 
for an isotropic work hardening material (2.13) reduces to [1] 


A (p , J, k) = 0. 


(4.1) 


To obtain the yield curve (3.10) in terms of- stress in- 
variants we use (3.6) ‘to arrive at 


where 


_P = 3> 2 

Po $ 2 +K£ 2 


(4.2) 


$ = /2/3 M , ? 2 = . 


(4.3) 


We see that (4.2) is of the form (4.1) with the strain hardening 
parameter -k identified with po . We can also write (4.2) in the 
form 


to = f(p; po) 

where To is the octahedral shear stress and 

f(p; Po) = "ti p(~ - i ) 1 / 2 
v^K p 


(4.4) 


(4.5) 


We also note that by using (2.16) and (3.12) we may write (4.2) 
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in the form 


6 - 


*-in (EisliSil, 


(l-n Q ) Y$2 


(4.6) 


It can easily be shown that the flow vector associated 
with the yield surface (4.2) is given by 


6 


'mr 


3 (1-n )p($ 2 -K£ 2 ) 


I • 

(P9(«)«„r +6K « m r> < E + — > <4- 


p -$2+K^2 


7) 


where 


gU) = <i> 2 - K(£ 2 + 6) . 


(4.8) 


Let us assume that there is no recoverable energy during 
shear distortion, then 


Y ’ - 1/2 0' 6 . 

’mr ' mr 


Also combining (2.24) and (2.26) we obtain 

0 ' = ~ * In (£) . 

(l-n 0 ) R 

Therefore 


0 = 


- — [ In (El^l+Kg jjj + Kin (2-)] 


(l-n 0 ) 




R 


(4.9) 


(4.10) 


(4.11) 


and 


1 mr 


3p ($ 2 -K£ 2 ) 


1 (pg(e)6 mr . + 6K a ) 0 11 + 5.(5 (4.12) 

mr mr _ mr 


3(l-n fl ) P 


Here 


6 " = ^ — (£ + . 


(l-n 0 ) p $ 2 +K£ 2 


(4.13) 
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To the above equations we add the equilibrium equations 

a , . . + X- * 0 (4.14) 

where is the body force, and the strain rate-velocity relations 

V = 1/2(v m,r + v r , m > (4.15) 

* 

where v^. = u r are the components of the velocity vector. Equa- 
tions (4.6), (4.12), (4,14) and (4.15) form a system of sixteen 
equations for the sixteen unknowns Y mr , 9"/ and v r . 
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5. PLANE STRAIN 


In this section and for the duration of the paper Greek 
indices will have the range from 1 to 2 . We further consider the 
case in which the soil constant tc is negligible in comparison 
with Xo - Under this assumption the elastic response is negligible 
and we may set 



•Consequently 


(5.1) 


y - y" 

mr 


(5.2) 


and we may drop the double primes to identify the irreversible 
strains. 

A state of plane strain is characterized by the assumption 


u a " V* 1 ' x 2 ) 

u 3 = 0 


(5.3) 


It follows that for plane strain 


Y ' - 0 , y — 0 

ms itt3 


(5.4) 


and therefore 


0 = y r 0 = y • 

’aa r aa 


(5.5) 


Now in view of (5.1) the constitutive relation (4.12) 


reduces to 
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Y =: r (pg(£)6 + 6Ka )6" 

^ 3p(4> 2 -K£ 2 ) ^ M 


(5.6) 


Also considering (5.4) we obtain from (5.6) the pair of equations 


a 13 = 0 23 * 0 


(5.7) 


pg(£) + 6Ka 3 3 = 0. 


(5.8) 


Equation (5.8) is a quadratic equation in cr 33 which may be 
solved to obtain cr 33 as a function of cr a ^. However, since g r 
a 33 must also satisfy (4.2) we may combine (4.2) and (5.8) so as 
to obtain 


3 3 


- 3 P<fr 2 ~2(d» 2 “3K)g 
2 ($ 2 +6K) 


where 


a ~ o 

aa 


(5.9) 


(5.10) 


With 'a 33 given by (5.9) the mean pressure is given by 


6Ktn+po * 2 
p = 2- 

2 ($ 2 + 6K) 


(5.11) 


Using (5.9)' and (5.11) we reduce (4.2) to 


Cjl 2 + C 2 ct 2 + C 3 p g a = l/4 - p 0 $ 


(5.12) 


where 


Ci 

C, 


K ($ 2 + 6K) 

($ 2 -3K) ($2-3K+9K ? ) 

($ 2 +6K) 

(1-K) ( $ 2 -3K) 


C, = -3 <t> 2 (1 + 


') 


(5.13) 


$ 2 + 6K 
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I 2 “ ° 0 °a ’ 

2 a0 0a 

Equation (5.12) describes the yield surface in plane-strain, 
Now from .(5.12) we find 


( 1_n o> (p„ 2 $ 4 -2C 3 ap 0 ) 


^ 2Cl °a0 + ^ 2< “ 2Cr+C3Po ^a0^ (5.15) 


where we have used (2.16) , (3.12), and (5.1). We can then show 
that with the yield locus (5.14), the associated flow rule is 


y = ± (2C ia -+(2C 2 atC 3 p 0 )6 ) 6 


a0 (C„a+2C 3 p 0 ) 


(5.16) 


where 


C 4 = 2Cj + 4C„ . 


(5.17) 


To the above equations we add the strain rate-velocity and 
equilibrium relations 


Y afS = 1 / 2 (v a,S + W 


(5.18) 


a + X 2 ~ 0. 

a0' 0 


Equations (5.12), (5. 16),. and (5.18) form a system of nine equa- 
tions for the nine unknowns pa, a , Y a g/ and v a - 
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6 p COMPARISON WITH THE "CAMBRIDGE THEORY"' 

In this section we shall compare the three-dimensional 
theory of Section 4 and the plane strain theory of the previous 
section with the corresponding theories presented by Roscoe 
and Bur land [2] . 

In order to make the comparison we must point out that ■ 
Roscoe and Burland [ 2 ] use only principal components of the 
stress and strain tensor. Using principal components of stress 
the invariant j is given by 

J = 1 /l[io x -a z ) 2 + (cr 2 -a 3 ) 2 4- (cr 3 - ai ) 2 ] (6.1) 

and is related to the stress parameter r used in [2 3 by 

r - . (6.2) 

Therefore from (4.3) 2 , £ = r/p and for K = 1, (4.2) reduces to 

2— = ^ (6.3) 

Po 4> 2 + K 2 

which is the equation for the yield surface of the three-dimen- 
sional theory given in [2] . 

Now combining (4.10) and 4.13) and setting K = 1 we obtain 

6 = 1 — + (6.4) 

(l~n 0 ) $ 2 + £ 2 p 

which is the equation given in [2] for the volumetric strain 
increment. 

0 

We consider the invariant e of the increment of irrecoverable 
deviatoric strain given by 
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(e ) 2 = e » e" . 

mr rm 

i 

Since we have assumed = 0 we may write this as 

(e) = £ £ . (6. 

mr rm 

Using principal components of total strain, e is given by 

, £ ~ — T Hyi-Ya ) 2 + (Y 2 -Y 3 ) + ( Y 3 “Y i ) (6.' 

vY 

o 

Roscoe and Burland [2] use e as a strain-increment parameter. 
Now from (4.7) we can show for K = 1 


p(4> 2 ~£ 2 ) 


and from this equation we find 


($ 2 -£ 2 ) 


Substituting. (4.13) into (6.3) and setting K = 1 we find 

e « — ( — (£ + . ’(< 

(1-n^) 4> 2 -£ 2 p $ ? +£ 2 

Equation (6.9) is the equation given in [ 2 ] -for the strain- 
increment parameter e . 


We observe that we may combine (6.7) and (6.8) so as to 


obtain 


E mr T mr r * 


Since e 1 = 0 it follows that 

rnir 


Y = x - + 1/366 . 

nir mr mr 
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Using principal components of stress and. strain (6.11) is equi- 
valent to the three equations 

• 

Y i = 1/3 [ (2a 1 -o 2 *-a 3 ) ~ + 03 

r 

• 

yV 1/3 [ tioa-os-oj) 0] (6.12) 

r 

• 

Y 3 ~ 1/3 I (2a 3 -a i -a 2 ) — + 0*3 

r 

which are the equations given in [2] for the three principal 
strain-increments. Thus, for K = 1 , the three dimensional theory- 
given in Section 4 is in complete agreement with the three-dimen- 
sional theory presented by Roscoe and Burland [2] for "wet" clay. 

We shall now compare the plane strain theory of the previous 
•section with that given in (2] . The plane strain theory given in 
[2] is derived under the assumption < = 0. This assumption was 
also made in the previous section. For K = 1, (5.9) reduces to 


= 3po $ 2 ~ 2 ($ 2 -3) g 

2 ($ 2 + 6 ) 


(6.13) 


which is the equation given in [2] for the determination, of arj 3 . 
Also for K = 1, (5.12) reduces to 

(4> 2 +6)X 2 + ($ 2 -3)a 2 - 3$ 2 p Q g - 1/4 p o 2 4> lt . (6.14) 

Using principal components of stress we reduce (6.14) into the 
'form 


(24> 2 +3) ( $ , 2 +a 2 ) - 3$ 2 p (a ,+a 2 ) + 2($ 2 -3)aia, = 1/4 p 0 2 4> 4 

12 (.6.15) 

which is the yield curve for plane strain in o 1 ~o z space given 

in [2] . 
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Roscoe and Bur land [2] introduce stress parameters t and x 
defined by. 

t = 1/2 (a +a 2 ) , t = l/2(aj-<7 2 ). (6.16) 

Using experimental observations as a guide Roscoe and Burla'nd 
[2J introduce the following approximation to- (6.15) > 


where 


,Po „1/ 2 

Oi - n(“- - 1) 

t 


h j = — , Q — -it = — — M. 

t /2 , /3 

The incremental stress-strain relations 
giyen in [2] are . 

s'. t 

(1-n ) fl 2 +co 2 t 


(6.17) 


(6.18) 

based on (6.17) , 

■ (6.19)- 


y = ... e 

(n 2 -w 2 ) 

where 

Y ~ Yj-Yz • (6.20) 

In order to compare our results with those in [2] we find 
it convenient to introduce Mohr's circle variables t, x, and xfs 
through 


t *= 1/2 (an + a 2 2 ) 

t * [1/4(01! - o 2 2 ) 2 + a 12 c 12 ] 1/2 (6.21) 

20 1 z 

tan 2\p =■ — — . 

G i 1 “ a 2 2 

In terms of stress parameters t and x the yield locus given by 



89 


(5.12) reduces to 

2Cj t 2 + 2 (C^+2C^) t 2 + 2C 3 p Q t = 1/4 p Q z $ 4 (6.22) 

where Ch (i=.. 1,2,3) are given by (5.13). For K = 1, (6.22) 
reduces to 

64> 2 1 2 + 2 ($ 2 +6) t 2 - 6$ 2 p 0 1 = 1/4 p,, 2 ^ ‘ (6.23) 

which is the equation for the yield locus given in- (2] in terms 
of stress parameters t and t. 

Equation .( 6 . 23) may be put into the form 

(1+Rl) _ El + (l - £ 2 J?.a . 2 > 0. (6.24) 

3 n 2 t 12t 2 


The expression (6.17) emerges from (6.24) under the following two 
conditions: (i) Q 2 /3 << 1 and (ii) ft 2 p 0 z /12t 2 << 1. It can be 

shown that under these conditions 


p<$ 2 - £ 2 ) ~ 2t(n 2 - w 2 ) . . (6.25) 

In view of (6.25) we obtain from (5.6) and (5.8) with K - 1 , 
the expression 


Y - - (a n - a 3 ,6 )0 

t (£J 2 ~u 2 ) aB 

Using this expression we construct 


<Y* - Y*2> - 8. 

:1 t(fi 2 -w 2 ) 


(6.26) 


(6.27) 


If we use principal components of stress and strain rate (6.27) 


will reduce to 




Equations (6.28) and (6.30) are in complete agreement with (6.19). 
Thus for K=l, the approximate plane strain theory extracted from 
the plane strain theory presented in the previous section is in 
complete agreement with that of Roscoe and Bur land [2] . 
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Figure 1. Triaxial Compression Yield Curve 
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SECTION IV 

AXISYMMETRIC PLASTIC FLOW 
1. INTRODUCTION 

The present investigation in theoretical soil plasticity is concerned 
with ideal soils whose postulated mechanical behavior is an approxi- 
mation to that of a wide class of natural soils. It is the object of 
this investigation to provide a theoretical analysis, valid under cer- 
tain mathematical and physical assumptions, that has applications to a 
fairly wide class of problems that concerns the general situation of 
quasi-static axially symmetric plastic flow. 


2. PRELIMINARIES 

We consider the particles of a dry soil mass continuum to be referred 
to a fixed rectangular Cartesian system x t (i = 1,2,3) and let u denote 
the displacement vector field. The components of the linear strain 
tensor Yi ^ are then given by 

2Yij = d jUl + S lUj (-2.1) 

We also let n, ho denote the current and initial porosity of the soil 
mass continuum respectively and define the change in porosity 5 by 

5 = n - no (2.2) 

which we assume to be of the same order of magnitude as the linear strain 


tensor Yi j . 
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In an earlier work we defined yield in a soil as a permanent irre- 
coverable deformation [l]. The strain may then be written as the sum of 
an elastic or recoverable part y* and a plastic or irrecoverable part 

YiV 

Yij » Via + yfi . (2.3) 

We shall also assume porosity to be given as the sum of an elastic part 
n 8 and a plastic part n p , 

n = n° + n p . (2.4) 

In view of this assumed resolution we can write 

g = g 8 + n p (2.5) 


where g 8 = n° - n 0 . 

The theory presented in [l] is based upon the hypothesis that the 
irrecoverable deformation of the soil is described by the pair (yf , n p ). 
Also in [l] we hypothesized the existence of a unique function P of 
plastic porosity n p and Terzaghi’s effective stress such that 

F(a u , n p ) - 0. (2.6) 

In plastic porosity-stress space (2.6) defines a six-dimens ionai hyper- 
surface 2 which was called in [ 1 ] , the state boundary surface. A soil 
particle is then said to be in a plastic state if the value of the stress 
and plastic porosity at the particle are such that (2.6) is satisfied. 

In the theory of soil plasticity presented in [l] a yield surface 
for the soil was defined to be the projection on stress space of curves 
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on £ along which the plastic porosity has the constant value n p . This 
is given by 


AOij, *) = 0 


(2.7) 


where 


H = H(n p ) 


( 2 . 8 ) 


is the strain hardening parameter. For a properly chosen yield sur- 
face we can associate a flow rule by invoking Drucker's postulates [l,2]. 


We assume that soil Irrecoverable deformations are mainly due to 
the relative motion of the soil solid particles and that during this 
motion the soil solid particles are essentially incompressible. It can 
then be shown that under this assumption 


0 p = - 


d-no) 


(2.9) 


where 0 P = y£ B 'is the plastic volumetric strain. .We note that in view 
of (2.4) it follows from (2.8) 

* = h( 0 p ), (2.10) 


i.e., the strain hardening parameter is a function of the irrecoverable 
volumetric strains. 

Until further notice we adopt the convention to consider compres- 
sive stresses as positive. We further consider the stress field for 
which 


a Br = 0, m # r 


( 2 . 11 ) 





95 


Such a' stress field is called a triaxial stress field. In a triaxial 
stress field it is convenient to use the set of generalized stress 
parameters (q,p) where 


q — d 3 a - 

p - '|'( a S3 + 2^11). 


( 2 . 12 ) 


The soil mass that we consider is one that obeys the Coulomb 
theory of internal friction according to which the strength of the soil 
is limited by its ability to resist shearing stresses. In the two 
dimensional q-p space the Coulomb failure surface reduces to two straight 
lines Fi and F 2 (Figure 1) meeting at a common point on the p-axis 
and whose equations are 


q = M(p + c cot<J>) (Fi) 
q = -M(p + c cot^>) (F 2 ) 


(2.13) 


where c and j b are respectively the unit cohesion and angle of internal 
friction and 


M = 


6 sino 
3 - sino 


1 - 6 sino 

3 + sino 


(2.14) 


The lines Fi and F 3 (Fig. 1) are respectively the compression and tension 
Coulomb failure lines. We introduce the parameter N defined by 


N = N(p) = M(1 cot<J>) (2.15) 

We shall use this parameter below. 

In [l] a one-parameter family of yield curves was constructed in q-p 
space for a cohesionless soil (c=0) under the assumption that M = M 1 . 
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We carry this assumption into the present study and consider the follow- 
ing curve in q-p space (Fig. 1), 


JE. a 
po 


Nu 3 

Nu 3 + KTf 


where T] = q/p and 


(2.16) 


N u = N(P U ) 

K = PQ_" . Pu 
Pu 


(2.17) 


Here p u is the value of the mean pressure at failure at a change in 
porosity § and po represents the pressure corresponding to § on the 
virgin isotropic compression curve of the soil. Equation (2.16) defines 
a one-parameter family of curves. For the parameter of this family we 
take po which we assume to be given by 

Po = Y exp (- ) (2.18) 


where Y is the initial yield pressure under isotropic compression and 3 
is a soil constant. We can also write (2.18) in the form 


po = Y exp (\0 P ) 


(2.19) 


where X = (l-no)/P. 

The yield locus given by (2.16) is in terms of triaxial stress field 
parameters (q,p) and as it stands it is only good for analysis under this 
particular stress field. Let us introduce the invariant 


J = T rfc Tj. u (2.20) 

where T rn is the stress deviator. We also introduce the parameter § 


through 
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P3 _ J 

^ ' P 3 


where now 3p = a aa . We can then show that the expression 


£ = 


$2 


po + Kg 3 


will reduce to (2.16) in a triaxial stress field. Here 


( 2 . 21 ) 


( 2 . 22 ) 


$ = /Z N u 
3 

Combining (2.22) and (2.19) we obtain 

6 p - _J | n [ £.(* 3 . + ] 

^ " (l-n 0 ) ln L 


(2.23) 


(2.24) 


The requirement that in stress space the strain-rate vector be 
normal to the yield surface leads to the following 


£ r a T%3 ..eS's" 1 r a ® P > 


' rQ p ($ 3 - kS s ) 


(2.25) 


where e|? a is the plastic strain deviator. Also from (2.24) we obtain 


qp - 3 r . P -1 

(1-no) *■ ^ - K§ a + P J 

where we have assumed K to be constant. 


2K§§ 


(2.26) 


We consider the stress states for which p = p u . We call such 
states, limiting stress states. For Oj j a limiting stress state (2.22) 
reduces to 


J ss p 3 $ 3 = — M 2 (p 4- c cotij>) 3 . 


(2.27) 


Also in this case (2.25), (2.26) reduce to 



98 


G P 

tif m 


_g.K_£ T 

(I-K) J r 


0 p 


0P = 


_J 

(1-no.) 


, 2K I 
Hl+K) 5 



(2.28) 


We recall that the theory of earth pressure is based upon the 
concept of states of limiting equilibrium satisfying Coulomb's law of 
failure [33. The defect of the theory, however, lies in its develop- 
ment without constitutive relations. Thus in the- theory a stress field 
can be found, in principle, without explicit knowledge of an acceptable 
velocity field. The need for the necessity of a compatible velocity 
field with a limiting stress field satisfying Coulomb’s law of failure 
led to the suggestion made by Drucker and Prager [ 4 ] of using Coulomb's 
failure criteria as a yield criteria and to treat the soil mass as a 
perfectly plastic material. Prediction of volume changes under the 
idealization of a flow vector normal to the Coulomb failure surface, 
however, were higher than those found by experiments. Below we shall 
state our concept of associating a flow rule compatible with limiting^ 
stress fields. 


Even though Shield [ 5 ] has shown that the interpretation of the 
Coulomb law leads to only one failure surface for three-dimensional 
stress fields we find that (2.27) gives a convenient valid generaliza- 
tion of the Coulomb law to three dimensions. Equation (2.27) defines 
■a surface in stress space that we call the limiting surface. The set 
of stress points that lie on the limiting surface are not on one yield 
surface. However, each limiting stress point does -lie on the curve 
defined by the intersection of some yield surface with the limiting - 
surface. Hence it is correct to associate a flow vector with each of 
the limiting stress points on the limiting surface. Here is the main 



99 


difference between our use of the Coulomb surface as a limiting surface 
rather than a yield surface. The strain-rate is then not normal to the 
limiting surface. 


3. EQUATIONS FOR AXIALLY SYMMETRIC DEFORMATIONS 

In a cylindrical coordinate system (r,9,z) we denote by (cr r , cr^, 

♦ * « • 

Tqj , T rzJ t P q) the components of the stress tensor, (e r , e^, e z , Ygz , 
• * 

Yr z j Yr q) the components of the strain-rate tensor and (u,v,w) the 

components of the velocity. We are interested in axially symmetric 

deformations in which the z-axis is taken as the axis of symmetry. 

Under the assumption of axial symmetry the shear stresses t qz , Lqj 

* • 

shear strain-rates Yq Z > Yr q and circumferential velocity v all vanish 
identically and the remaining stresses, strain-rate components and 
velocities are only functions of (r,z,t). 


The strai'n-rate-velocity relations for axial symmetry are 


• du 


* u 

e 0 " r 

» dw 


• l ,3u , dw. 

Yr* “ 2 + 


Also for quasi-static conditions the stress components satisfy the 
equilibrium equations 





= 0 


dT 

■5? 


dcr, 

dS 


rz +^+ - rz + P g 


(3.2) 



PRECEDING PAGE BLANK NOT FILMED 
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assume the soil mass to be in a limiting state of equilibrium so that in 
effect the stress field must satisfy (2.27) in addition to (3.2). 

In view of (3.6) we may reduce (2.27) into 
s = f (a) = a a sin <j> + a c costf 

or ' (3;8) 

s “ g(p) = b p sin 4 + b c coscj> 


where 


_ 3 _ a 

a 3 - (1 + 3a) sintj* ’ 1 + Ct a sincj) 


(3.9) 


We can write cx r , cr z , T rz in terms of the two independent variables s 
and f 

cr r = — (s csc<j> - ac cot<J>) + s cos2i}r 

3 

- a z = — (s csc4 - ac„cot(j>) - s c os2ty (3.10) 

3 

T rz <= s sin2'li. 

Also we can show that under the Haar-von Karman hypothesis 


- Qq = s(cos2f - 3or) 


(3.11) 


Substitution of (3.10), (3.11) into (3.2) leads to the two differential 
equations for the two unknowns s and 

(cscijl + a cos2t) ■§” + a sin2ty 

- 2as C sin2<|i - cos2t(r (coa2t|» - 3tf) ] = 0 


a sin2ty + (csc4 - a cosZ^) 


(3.12) 



102 


+ 2as [ cos2'|' + sin2<(i ~ sin2i|» ] + pag - 0 


Now from (2.28), (3.6), and (3.10), we obtain 


'<•’ ~ 4(1 3 - K) f 811,2111 ® 


6 = -j-, [2K - + (1 

(1-no) (1+K) s v 


(3.13) 


K) £ ] 
-P 


where we have dropped the identification of the plastic strains since 
we are neglecting elastic strains. 

The volumetric strain-rate 9 is given by 


a _ Sw u 

e --^ + ^ + 7 


(3.14) 


Combining (3.13), (3.14), and (3.1)4 we obtain 


. bn bu . 3w r u 

sin2ty igp - h, + sin -g— + sin2'lf — = 0 


(3.15) 


where 


» • Ms) - 2<1 - K & — 


be cosqi) 


(3.16) 


Equations (3.12), (3.15) together with the equation of istropy (3.3) 2 


6u 3u Ow 


dw 


tan2 * £-■£ -£■- ^5=0 


(3.17) 


form a system of four equations for the four unknowns s,i)J,u, and w. 
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4. STRESS AND VELOCITY CHARACTERISTICS 


We examine the case for which 3 - -1, i.e., cr 3 = cr 2 . For this case. 


a = 1 

k 3 - sino 


(4-1) 


Using known methods we may obtain the equations for the character- 
istics of the system of equations governing s and tjf. The slopes of 
the characteristics in the z-r plane are given by 

dz. _ sin 2 ty ± cose ,, 

Mr 1,3 cosgty + sino. 


We see that the characteristics will be real and non-zero. It follows 
that the system will be hyperbolic. It is convenient to introduce the 
angle 9 through 



(4.3) 


and to name the characteristic with slope tan(l|i - cp) an ^-line, and that 
with slope tan(ij) + 9 ) a 3 - line, thus 


If or = tan ^ " ?) 

P = tan(ijJ + cp). 


(4.4) 


The a and 0 characteristic directions are illustrated graphically in 
Fig. 3. 


The equations along the characteristics can be shown to be 
1 2 s 

cotqids - 2 sdty - (pg cos((f+cp) - — sincp cosi|f)dS a = 0 on cr-line 
1 2 s 

cotq>ds + 2 sdf + (pg cos(i)J-cp) + — sincp cos^dS^ = 0 on 3 -line 


(4.5) 



104 


where S- and S Q are arc lengths measured along the a and 0 lines 

Of p 

respectively. Introducing the quantity X defined by 


X “ cotdi In — 
So 


(4.6) 


where sq is some reference stress, we may write (4.5) in the form 


dX - 2di|r + ^ (cost)) dr - (1 - sincj>)dz) 

- (sin(|> dr - cos(j> dz) exp (-X tancj>) - 0 an tf-line 

So 

(4.7) 

*• 

dX + 2di|r + — (cos|> dr - (1 - sinc^)dz) 

+ (sin<[ dr + cos^ dz) exp (-X tan(j>) = 0 on 0-line. 

So 

We consider the case when if has a constant value if along a stress 
^-characteristic (say). In this case (4.7), becomes 

- B exp (-X tan<j>) + ^ = 0 (4.8) 

where 


A = cos(ji - (1 - sintj>) tan(if o - <p) 

B = (siniji - cos(|p tan(i|l 0 - cp)) — 

So 

The general solution of (4.8) can be shown to be 

__ KS 0 ■ Bs 0 tan|> 

~ r A tanq> (1 + A tan<ji) 


(4.9) 


(4.10) 


where K is a constant of integration. At a later time we shall refer 
to (4.10). 


The slopes of the characteristics associated with the system (3.15), 


(3.17) are given by 
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dz 

dr 


1 13 


sin 2 i|J ±/l-K, g 
cos 2 t[f - K 


C4.ll) 


where now in view of (4.1 ) 3 


* - ^ ■*«* S 


(4.12) 


We assrime h < 1. This being the case the system (3.15), (3.17) will 
be hyperbolic with characteristic directions given by (4.11). By 
defining the angle T] through 


- sin 2 T 1 = h 


(4.1-3) 


we may write (4.11) in the form 


(f“)i “ tan(t - 9) 
(•3^)2 53 tan(t + 9) 


(4.14) 


where 


^ TT „ 

9 = T ' ^ 


(4.15) 


Comparing (4.15) with (4.3) we see that in general the velocity charac- 
teristics do not coincide with the stresB characteristics. The equations 
along 'the velocity characteristics are 




dz. 


- (cos2ty - sin2'n) _i (sin2ijf - (—) 1 , a cos 2$) - dz = 0 


dr 


r 

(4.16) 
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5. INDENTATION OF SEMI- INFINITE SOIL MASS BY A 
LUBRICATED CIRCULAR RIGID CONE 


As an application of the theory presented we consider a problem 
that is related to the cone- penetrometer test used in soil mechanics. 
This problem is associated with the incipient plastic flow in a semi- 
infinite region of soil due to load applied through a lubricated circu- 
lar rigid cone. 

We define the origin of cylindrical coordinates as shown in Fig. .4 
and we shall suppose the soil to occupy the semi-infinite region 
z ^ 0. In addition to the load applied to the soil through the circular 
cone we take into account the normal stress t (Fig. 4) which may cor- 
respond to atmosphere pressure or an equivalent surcharge. 

The boundary conditions on f, s, and velocity components for the 
problem of interest are 


w = 
u = 
* = 

s = 


6 

const. 

const. 

0 


on z = (R-r) cot§ 


(5.1) 


t sine -f c cosn 
1 - sino 


on r = R, z = 0 


(5.2) 


The boundary conditions (5.1), (5.2) together with the governing equa- 
tions (3.12), (3.15), and (3.17) define a boundary value problem. We 
note that the velocity boundary conditions are not sufficient to com- 
pletely determine the velocity field anywhere in the soil mass. This 
indicates that we may impose further arbitrary conditions on the solu- 
tion to find an acceptable velocity field. However we should not expect 



107 . 


this solution to be unique. Cox, Eason, and Hopkins [6] discuss the 
uniqueness of a similar problem. 

Following Cox, Eason, and Hopkins [6] we expect the stress charac- 
teristic net to exhibit the geometrical features depicted schematically 
in Fig. 5. The unknown of interest is the limiting mean pressure at 
the face of the cone. To find this we need only to consider that part 
of the stress field bounded by 0A, AB and BCDO. Lines such as PS will 
be cf-lines while those such as PQR will be {3-lines. In particular 
ODCB is the 3- line through the apex of the cone 0. We shall discuss 
this 3- line below. 

A solution for’ the stress field can be constructed following argu- 
ments similar to those used by Cox, Eason, and Hopkins [6]. Thus from 
a knowledge of s and f on AB, s and can be determined on ABC by 
using (4.4) and (4.5)-. At r = R, z = 0 a singularity is introduced 
at which s and will be multivalued. This fact together with the known 
values of s and ijf can be used to determine s and tjf on ACD. Finally 
the now known values of s and f on AD together with the known values of 
'if on 0A determine s and l|r on ADO. From this solution the limiting 
stresses at the face of the cone can be determined. If we find a velo- 
city field compatible with this stress field then the solution will 
given an upper bound for the limiting stress field. If this stress 
field can be extended in such a way so as to satisfy the conditions 
of limiting equilibrium then the solution is said to be complete and 
will be a lower bound for the limiting stress field, [6]. 

Now by considering the singular point A as a 3-line of zero 
length we, find the conditions at point A, 
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0 £ iji <; 6 



(5.3) 


In order to extend the solution into the rigid region we must 
know the location of the boundary between the rigid and deforming regions. 
Since at this boundary the velocity components or their spatial deriva- 
tives must have some discontinuities it follows that the boundary between 
rigid and deforming regions must be a velocity characteristic line. 
Considering that in the theory presented here the velocity and stress 
characteristics do not coincide we reason that the stress P-charac- 
teristic curve ACDO is not the boundary between .the rigid and deforming 
region. 


Let us name the first velocity characteristic an o' -line and the 
second velocity characteristic a P'-line. In Fig. 4, the curve B'O 
depicts the velocity 3* -line through the apex of the cone. Here we 
have assumed that 



(5.4) 


Now the velocity field must accomodate the incipient motion of the 
rigid cone. For this reason AO must lie within the deforming region. 

The simplest configuration that can occur is when the P'-line B'O is the 
boundary between the deforming and rigid region. We shall return to 
this question below. 

The stress field and stress characteristic net in the region BAOD 
can be determined by numerically integrating (4.4) and (4.7). Using 
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pgR as the reference stress s 0 indicated in (4.6) we replace (5.2) and 
(5.3) by 


v|f * 0 

i | .t* sin® + c* cos® 

x ' ct "* |n < — r- ,14 — 1 

i’ ) ,t* sin® + c* cos 

* = cot, » ln < — rrife| 

0 ^ \|( ^ 6 

Here t* = t/pgR, c* - c/pgR. The numerical integration is based upon 
the approximation of (4.4) and (4.7) by finite difference equations [7]. 

It remains to determine the velocity of 3 '-line B'O which we have 
assumed to be the boundary between the deforming and rigid region. Since 
in addition to i|f, X is now known on. AO we can determine the slope of 
B'O. at 0. Suppose we extend the stress a- lines beyond BCDO as straight 
lines with slopes obtained by using the calculated values of i|i on 
BCDO. Then along each of these stress Qf-lines (4.10) holds with the ar- 
bitrary constant K evaluated for each Of- line by values on BCDO. With 
reference to Figure 6 we can determine the coordinates (r T , zj) of • 
the point T of the intersection of the straight stress o'- line with 
the velocity P'-line as follows: We assume values of the slopes known 

at S and M. Then the first approximation (ri, zi) to (r T , z T ) is 


on r > R, z = 0 
, } (5.5) 

® at r = R, z = 0 


obtained from 


zi - z fl = (ri - r B ) tan( - «P) 
z x - z„ = (n - r 3 ) tan i|r* H 


(5.6) 
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where is known with the calculated value. With this calculated 

value of r* we obtain through (4.10), (4.12) the first approximation 
Hi to h t . With this value h, we can compute the first approximation 
t|(*i to iJf* T . We can then use the average value of f*! and i(i* M to obtain 
second approximations. Using this iteration procedure we can deter- 
mine the location of the p'-line B'O and the value of S at each point of 
intersection of the straight stress a-iine with the P'-line B'O. 

Once the stress distribution along B'O is known we can determine 
the total vertical load P, exerted by the indenter, by simple statics. 
The ultimate bearing capacity qu is then obtained from 



6. CONCLUSIONS 

We -have presented an approach of obtaining ultimate loads that 
differs from the traditional approach used in soil mechanics. The tradi- 
tional approach for obtaining ultimate loads for soil mechanics problems 
is to exhibit an equilibrated limiting stress field, solve for the load 
equilibrating the stress field, and- then simply term this load a 
"failure" or "ultimate" load. The approach presented here, however, 
demands more of a load before it is termed an ultimate load. The ulti- 
mate load must, in addition to equilibrating a limiting stress field, 
be associated with a deforming solution in a theory involving material 
deformation. 

The problem formulated in Section 5 can be tied to the cone pene- 
trometer test. The amount that the cone has penetrated is equal to the 
amount h s that the shaft has penetrated plus the height h 0 = R tan6 
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of the cone.- The distance h a is reflected in the equivalent surcharge 
t = pgh e which is a boundary value. We can then associate with an 
ultimate load P a depth h = h c + h e . Thus we can plot a cone load 
P vs. depth h curve. This, of course, is for a given soil with cohe- 
sion C, angle of internal friction <j> and bulk mass density p. 

We introduce the quantities q u *» h*, and c* defined by 


qu* 


= HlL. 

pgR 


h* 



c* 


c 

PgR 


( 6 . 1 ) 


Then for a ■ given soil with cohesion c, angle of internal friction tj>, 
bulk mass density p and given fixed cone geometry R and 6 we can find 
an ultimate bearing capacity q u at depth h. We can then generate the 
curve 


q u * = f (h*; c*, <j) (6.2) 

wherein | and c* appear as parameters. 

The solutions for q u * most widely used in soil mechanics are 
of the form 


q u * = c*H c + h*H q (6.3) 

where H c , H q are bearing capacity factors which depend on (j> alone. It 
is possible to obtain from (6.2) values of N c , N q such that the form 
(6.3) holds. However by doing so we would have to accept the obvious 
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consequences of the form (6.3), viz., linearity in c* and h*. 

The use of the cone penetrometer as a useful test of obtaining 
in situ properties of soils is highly dependent on being able to solve 
the following problem: given the response curve q„* vs. h* determine 

the parameters c* and that correspond to this* curve. In this section 
•we have presented a theory to generate q u * vs. h* curves for different 
values of c* and <j>. It remains to look into a systematic and logical 
manner of using response q u * vs. h* curves, to obtain values of c* and 

<j). 
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Figure 1. Triaxial Compression Yield Curve 



Figure 2. Directions of Principal Stresses 
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